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ABSTRACT 

We present MIPS 24/im and 70/im photometry for 205 members of the Upper Scorpius OB As- 
sociation. These data are combined with published MIPS photometry for 15 additional association 
members to assess the frequency of circumstellar disks around 5 Myr old stars with spectral types 
between BO and M5. Twelve stars have a detectable 70/im excess, each of which also has a detectable 
24/im excess. A total of 54 stars are identified with a 24/im excess more than 32% above the stellar 
photosphere. The MIPS observations reveal 19 excess sources - 8 A/F/G stars and 11 K/M stars 
- that were not previously identified with an 8/tm or 16/im excess. The lack of short-wavelength 
emission and the weak 24/im excess suggests that these sources are debris systems or the remnants of 
optically thick primordial disks with inner holes. Despite the wide range of luminosities of the stars 
hosting apparent debris systems, the excess characteristics are consistent with all stars having dust 
at similar orbital radii after factoring in variations in the radiation blowout particle size with spectral 
type. The results for Upper Sco are compared to similar photometric surveys from the literature to 
re-evaluate the evolution of debris emission. After considering the completeness limits of published 
surveys and the effects of stellar evolution on the debris luminosity, we find that the magnitude of 
the 24/im excess around F-type stars increases between ages of 5 and 17 Myr as found by previous 
studies, but at <2.6<r confidence. For B7-A9 and G0-K5 stars, any variations in the observed 24/im 
excess emission over this age range are significant at less than 2a confidence. 

Subject headings: open clusters and associations: individual(Upper Scorpius) — planetary sys- 
tems:protoplanetary disks, debris — stars:pre-main sequence 



1. INTRODUCTION 

Dusty debris observed around main sequence stars is 
a potential diagnostic of planetary systems. The cur- 
rent paradigm is that the debris originates from the col- 
lisional grinding of planetesimals to micron-sized grains, 
where the collision rate is dictated by gravitational 
interactions between planet s and a planetesimal belt 
(Williams & Wctherill 1994). If this hypothesis is cor- 
rect, the spatial asymmetries in the debris that have been 
revea led by high resolution imaging (e.g. iHolland et alj 
1998) encodes information on the eccen tricity, mass, 
and even migration of the orbiting planets (jLiou fe Zookl 
119991 iWvattj 12003). Unfortunately the resolution and 
sensitivity required to image the debris can be achieved 
for only a few disks with current instruments. 

A complementary approach to detailed studies of in- 
dividual disks is to study the ensemble properties of de- 
bris systems. Debris disks have bee n studied extensively 
over t he entire sky u s ing I R AS (iBackman fe Paresca 
[1991 ILagrange etall l2000t IRhee et all 120071) and 
in targeted regions using ISO (lHabing et all 120011 : 
iSpangler et"aT1l200l IDominik fe Decinl 120031) The 
Spitzer Space Telescope ( Werner et al.l [2004D has ex- 
panded on these studies to produce a vast database of 
debris disks (e . g. iRieke et al.||200"5 jBrvden et al J 12006 ; 
Su et al J 12009: IGautier et al.M2007t ICurrie et all 12008a 



Carpenter et al J 120091 see also reviews by iWerner et al.l 
20061 iMever et al.l 120071 and iWvattl |2008|) that ~erT 
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compass a broad range of spectral types (B-M stars), 
environments (clusters, associations, field stars), and 
ages (3 Myr to 10 Gyr). Comparison of these 
data to theoretical models have yielded insights on 
the planetesimal be lts that produce the debris dust 
(jWvatt et al.ll2007bf). the collisiona l history of planetes- 
imal belts (IDominik fe Decinl 120031 1 Wvatt et alj|2007at 
lLohne et al. 2008[) and the for mation of planetary sys- 
tems ( Kenvon fc Broml ey 2008, and references therein). 

A critical issue emerging from Spitzer studies is es- 
tablishing when the debris phenomenon is initiated. In 
the Kenyon & Bromley models, a planetesimal belt pro- 
duces low levels of debris emission in the early stage of 
planetary accretion. At a given orbital radius, dust pro- 
duction reaches a maximum when 1000-3000 km sized 
bodies form and ignite the destructive collisional cascade 
(|Kenvon fc Bromlevl 120041 ). As planet formation propa- 
gates through the disk, the disk is depleted of planetesi- 
mals through repeated collisions and the debris produc- 
tion eventually declines. In theory, the evolution of the 
debris emission can constrain the formation time of plan- 
ets and the c ollisional cascade. Indeed , bas ed on Spitzer 
observa tions, IHernandez et al] (|2006f ) and ICurrie et al.l 
(2008aJ) have suggested that the debris luminosity peaks 
around 10- 30 Myr for A-F st ars before declining toward 
older ages (jRieke et alJl2005f ). 

Additional observations are needed to validate 
previous conclusions since few surveys for debris 
emissio n at ages of < 10 Myr have been com- 
pleted (IHernandez et all l2006t iCarpenter et ail 120061 : 
IGautier et al.ll2008t ICurrie et all l2009aT) . The Upper 
Scorpius OB Association (hereafter, Upper Sco) is an 
important region in studies of disk evolution since it 
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provides a snapshot of disk properties at an age of 
~ 5 Myr when most opt ically thick disks have dissipated 
(Hernand ez et al] [2007) . Spitzer observations can then 
probe the onset of debris production as planets form. 
Since association members have been identified from B to 
M stars, disk properties can be investigated over nearly 
two orders of magnitude in stellar mass. Fin ally, at a 
mean distance of 145 pc (|de Zeeuw et al.lri999D . Upper 
Sco is more than a factor of two closer than other clus- 
ters or associations of similar age. These traits permit 
sensitive photometric studies of debris properties. 

In previ ous studies, we investig ated the 4.5-16/im pho- 
tometric (|Carpenter et al.l 120061 hereafter 'Pape r I) and 
5-35/mi spectroscopic (|Dahm fc Carpenterll2009l ) proper- 
ties of circumstellar disks in Upper Sco. The principle 
result from these studies is that the circumstellar disks 
detected at wavelengths shortward of 16/mi show a di- 
chotomy with spectral type. Circumstellar disks around 
late-type stars (K and M spectral types) are likely opti- 
cally thick primordial disks formed as a result of the star 
formation process, while early-type stars (B and A types) 
appear to be surrounded by debris disks. Surprisingly, 
no disks were detected around the F and G stars. 

The photometric data presented in Paper I spanned 
wavelengths between 4.5/im and 16/im. To probe for 
cooler dust, we present here 24/im and 70/im photomet- 
ric observations of 205 Upper Sco members. These new 
Spitzer observations of Upper Sco are presented in Sec- 
tion [21 In Section (3[ we use these data to identify stars 
that have excesses in the 24/im and 70/im bands. The 
nature of the excess sources are investigated in Section [J] 
and properties of debris systems around stars with a wide 
range in mass are compared in Section In Section [6[ 
we compare the Upper Sco results with published Spitzer 
surveys to investigate the evolution of debris systems. 

2. OBSERVATIONS AND DATA REDUCTION 

The selection criteria for the Upper Sco sample are 
described in Paper I. Briefly, the sample was con- 
structed from membership lists estab lished by Hippar- 
cos astrom etrv (Ide Zeeuw et al] 119991). color- magnitude 
diagr ams (jPreibisch fc Zinneckerl Il999t iPreibisch et al 
2002|), and X-ray surveys ([Walter et al.l 119941: iMartm 



1998; IPreibisch etHHlOM iKunkel 119991 : iKohler et al 



20001 ) . Since membership selection was based upon stel- 
lar properties (proper motions, photospheric colors, X- 
ray activity) unrelated to circumstellar disks, the sample 
is believed to be unbiased with respect to the presence or 
absence of disks. From the compiled list, we selected 205 
stars 4 to sample spectral types between B0 and M5. This 
sample was supplemented with 15 Upper Sco members 
with G/K spectral types that were observed by the FEPS 
(Formation and Evolution of Planetary Systems) Spitzer 



4 IBouv fc Martin (2009) found that the star [PBB2002] USco 
161021.5—194132 deviates from the mean proper motion for Up- 
per Sco and may be an interloper. The reduced x 2 from the 2MASS 
point-spread-function fitting photometry is the second highest 
among the nearest 100 stars to this source that have a 2MASS qual- 
ity flag of AAA. The astrometry may therefore be contaminated 
by a binary. The luminosity of this st ar is also higher than the 10 
nearest stars of the same spectral type (Prcibisch ct al. 2002), con- 
sistent with a binary contaminant. Given that the lithium equiv- 
alent width (Prcibisch ct al. 2002) and apparent brightness are 
consistent with Upper Sco membership, we retain the source in 
our sample. 



Legacy program (|Mever et al.ll2006[ ). The FEPS targets 
were selected based on similar criteria as our study and 
should also be an unbiased sample with respect to the 
presence or absence of disks. 

Other Upper Sco members have been observed with 
Spitzer but are not include d in our sample. The Her- 
big Ae/Be star HIP 79476 (jHernandez et all 120051 ) was 
observed with the Spitzer Infrared Spectrograph (IRS), 
but not with MIPS. We did not include the G6/G8 star 
HD 143006, which was also observed by the FEPS pro- 
gram as part of survey of gas in disks. HD 143006 was 
recogni zed as a Upper Sco member based on an IRAS 
excess (|Odenwaldl Il986l ) and thus would bias the sam- 
ple. We als o do n ot use the three B stars observed by 
Rieke et al.l (120051) and the five F/G stars observed by 
Chen et all ([2005) since IRAC data are not available for 
these sources, which is needed for the analysis presented 
in Section [3l Finally, we did no t include brown-dwarfs 
in Upper Sco (|Scholz et al.ll2007D . 

In Paper I, IRAC (4.5/im and 8/im) and IRS peak- 
up (16/im) photometry was presented for 218 of the 220 
stars. IRAC photometry for one additional Upper Sco 
member observed by FE PS (ScoPMS 52) is now available 
([Carpenter et al.ll2008f ). The remaining star (HIP 80112) 
was not observed with IRAC since it would have satu- 
rated the detector. The 4.5/im flux density for this source 
was estimated from groun d-based M-band observations 
([van der Bliek et al.lll99l . 

We present here MIPS 24/im and 70/im photometry 
for 205 stars in Upper Sco. The data reduction proce- 
du res follow that us e d by t he FEPS program as described 
in iCarpenter et al] (|2008l ) , and a only summary of the 
salient steps is provided here. Table [T] presents the 24/im 
and 70/im MIPS photometry for the 205 stars observed 
for this study. MIPS photometry for th e remaining 15 
stars are listed in ICarpenter et afl (|2008f ). 

2.1. MIPS 24fim 

MIPS 24/im observations for 205 stars were obtained in 
photometry mode. The exposure time per image (either 
3 or 10 sec) and the number of dithered images were set 
to achieve a signal-to-noise ratio of at least 10 on the 
expected photospheric brightness. MIPS 24/im images 
were processed starting with the Basic Calibrated Data 
(BCD) products produced by the Spitzer Science Center 
(SSC) pipeline version S17. Individual BCD images were 
combined to derive a flat field that removed long-term 
gain changes in the MIPS array. Flat field images were 
derived from stars that are not surrounded by nebulosity. 
If nebulosity is present, a flat-field image from another 
star was used that was observed close in time and with 
the same exposure time. 

Point (source)-response- function (PRF) photom- 
etry was performe d wit h the MOPEX package 
(jMakovoz k Marleaul [2005). The empirical PRF 
distributed with MOPEX was fitted to the individual 
BCD images simultaneously (as opposed to the mo- 
saicked image) using a fitting area of 21x21 pixels (1 
pixel « 2.55") for most sources. A 5x5 pixel fitting 
area was used for sources that have spatially variable 
nebulosity near the point source position. From visual 
inspection of the mosaicked images, the PRF of the 
Upper Sco target overlaps occasionally with that from a 
nearby source. These contaminating sources were fitted 
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with a PRF simultaneously. The measured 24/xm flux 
densities and internal uncertainties from the PRF fit 
are presented in Table [TJ The images were processed 
wit h a calibration fact o r of .0447 MJy sr _1 . Follow- 
ing lEngelbracht et aT] (|2007f ). we adopt a calibration 
uncertainty of 4%. 

Cirrus and extragalactic sources may contaminate the 
MIPS photometry and create the appearance of an in- 
frared excess. Since we expect the emission from a cir- 
cumstcllar disk to be point-like and centered on the star 
at the distance of Upper Sco, contaminated 24/im pho- 
tometry can be identified from emission that is extended 
or offset from the stellar position. 

The maximum astrometric offset observed between 
2MASS and the Spitzer 24/im position is 1.6" for 
HIP 78820. The dispersion in the offsets is the same 
(~0.4 to 0.55") for stars with and without infrared ex- 
cesses (see Section 13. I|) . indicating that any contam- 
ination is not systematically producing offsets in the 
measured coordinates. This empirical result is consis- 
tent with the cont amination anticipated from extragalac- 
tic source counts (|Papovich et al.l 120041 ). For each star 
in the sample, we computed the expected number of 
galaxies located within the full-width-at-half-maximum 
(FWHM) size (6") of the MIPS 24/im PRF that can pro- 
duce a 30 % ph otometric excess at 24/im. (As described 
in Section 13.11 30% is approximately the minimum ex- 
cess detectable with these observations.) We find that 
~ 1 star may be contaminated by an extragalactic source 
bright enough to produce such an excess. 

A more significant level of contamination may arise 
from cirrus and nearby stars. We used both quantitative 
techniques and visual inspection of the images to iden- 
tify such cases. First, we compared PRF photometry 
to aperture photometry obtained with a 6 pixel (I pixel 
= 2.55") aperture diameter, which is approximately 2.5 
times the FWHM size of the MIPS 24/rni PRF. For 38 
out of the 220 stars in the sample, the aperture photom- 
etry deviated from the PRF photometry by more than 
10%. We chose 10% as the threshold since such devi- 
ations would impact the significance the 24/im excesses 
identified in Section [3.1l Visual inspection of the images 
revealed that in the vast majority of cases, the differ- 
ence between the aperture and PRF photometry is most 
likely due to low signal-to-noise since the aperture size 
is much larger than the optimal extraction size for rela- 
tively faint sources (|Navlorl ll998). The photometry was 
deemed questionable for 6 stars based on the presence of 
extended nebulosity (HIP 77840, HIP 78265, HIP 79404, 
and [PBB2002] USco J161052.4-193734) or because the 
24/tm emission, while centered on the stellar position, is 
extended (HIP 77858 and HIP 80338). In addition, for 
the star ScoPMS 17, a known contaminating source is re- 
solved in the 16/im images (see Paper I) but is unresolved 
at 24/im. These seven stars are flagged in Tabic [TJ 

2.2. MIPS 70nm 

MIPS 70/im observations were obtained in photome- 
try mode with an exposure time of 10 sec and the small 
field size dither pattern. The number of cycles was fixed 
at 4 for all stars. MIPS 70/im images were processed 
with SSC pipeline version S17 that removes the bias, 
subtracts a dark image, applies a flat field correction, 
and linearizes the pixel response. Individual BCD images 



were mosaicked with the Germanium Reprocessing Tools 
(GeRT) software package S14.0 version 1.1 developed at 
the SSC. A 40" x 40" region centered on the source posi- 
tion was excluded when computing the column and time 
filtering such that the filtering process is not biased by 
the presence of a bright source. Filtered images were 
formed into mosaics with MOPEX (|Makovoz fc Marleaul 
120051) . 

Aperture photometry was performed on the 
MIPS 70/im mosaics with a custom version of IDLPHOT 
using an aperture radius of 16" (4 pixels on the coadded 
images). The sky-level was computed as the mean 
pixel value in a sky annulus that extended from 40" to 
60". The aperture was centered on the expected stellar 
position computed from the world coordinate system 
keywords contained in the FITS image headers. No 
centroiding was performed since the signal-to-noise ratio 
of most 70/xm measurements are less than 3. Visual 
inspection of the 70/im mosaics identified 19 images 
where a point source was located within the outer 
sky annulus or the aperture radius, but offset from 
the 2MASS stellar position by more than 4". Before 
measuring the aperture photometry, a PRF was fitted 
to the apparent contaminating source and subtracted 
from the image using MOPEX. 

The MIPS 70/im photometry and internal uncertain- 
ties are presented in Table [TJ Stars for which the 70/im 
photometry was measured on PRF-subtracted images are 
marked in the table. The adopted calibration factor is 
702.0 MJy sr" 1 / (DN s" 1 ) with an uncertainty of 7% as 
reported on the SSC MIPS calibration web pages. 5 

Of the 220 stars in the sample, 15 have a 70/im 
signal-to-noise ratio > 3. Photometry for two of these 
stars is compromised. ScoPMS 17 is contaminated by a 
nearby source resolved in the 16/im images (see Paper I). 
HIP 80338 is surrounded by nebulosity and it is ambigu- 
ous if the 70/im emission is associated with the star. In 
the remainder of this paper, we do not include these two 
stars in analysis of the 70/im data. 

3. SOURCES WITH INFRARED EXCESSES 

We now combine the 24/im and 70/im photometry for 
220 Upper Sco members (i.e. 205 stars from this study, 
and 15 st ars from FEPS; see S ection HD with 2MASS 
J and K a dSkrutskie et al.l|200ll and IRAC 4.5/im (Pa- 
per I. lCarpenter et al.ll2008 ) photometry to identify stars 
with infrared excesses. All ph otometry was deredde ned 
using the extinction law from iCardelli et al.l ( 1989) as- 
sumin g Ry — 3.1 for A < 2/im. the iChapman et al.l 
(2009) reddening law for I RAC and MIPS 24/im, and the 
reddening law compiled by Mathis ( 1990) for 70/im. The 
visual extinction toward individual stars were obtained 
from the literature, or derived by us from published spec- 
tral types and optical/2MASS p hotometry following th e 
general procedure described in lCarpenter et all {2008). 
The visual extinction is less than 2 mag for 97% of the 
stars, which implies a correction of less than 4% when 
dereddening the 24/im photometry. The adopted spec- 
tral types and visual extinction for each star are listed in 
Table El 

3.1. 24-[im excesses 

5 http:/ /ssc. spitzer. caltech.edu/mips/calib 
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Figure [JJ presents the 24/im to 4.5/im flux ratio (= 
^24/4.5) as a function of dereddened J — K s color for the 
Upper Sco sample. The J — K s color should represent 
the stellar photosphere for most sources since only two 
sources have an apparent -ftT s -band excess (Paper I). We 
chose to normalize the 24/im photometry to the 4.5/im 
band as a compromise between reducing the uncertain- 
ties in extinction corrections and the range of intrinsic 
photospheric colors (favoring longer wavelengths), and 
biasing the results due to weak emission from circum- 
stellar dust (favoring shorter wavelengths) . The locus of 
points with log -R24/4.5 ~ —1.4 have flux ratios roughly 
consistent with the stellar photosphere for emission in the 
Rayleigh- Jeans limits. In practice, the value of R24/4.5 
varies systematically with (J — K s ) , indicating changes 
in the photospheric color with spectral type. A number of 
sources have values of -R24/4.5 that lie substantially above 
this locus. These sources include the 35 stars from Pa- 
per I that have an 8/im and/or 16/im excess as indicated 
by the black circles and crosses. 

We adopted an iterative procedure to determine the 
locus of photospheric colors in the -R24/4.5 vs. J — K s 
diagram. The trend between log(i?24/4.5) and (J — K S ) Q 
color was fitted using robust linear regression after first 
excluding the known 35 stars with 8/im or 16/im excesses 
and the 7 stars flagged with suspect 24/im photometry. 
Sources that deviated from the best fit line by more than 
three times the dispersion of the residuals were removed, 
and the fit was repeated until no additional sources were 
excluded. 

The best-fit relation is shown as the solid line in Fig- 
ure [TJ The dispersion about the best-fit line is a = 
0.03 dex (i.e. 7% of the photosphere) for the sources 
included in the fit. The largest deviation below the best- 
fit line is — 2.5a. Given the sample size of 220 stars, we 
expected only one source more discrepant than — 2.5c 
for gaussian noise. The dispersion of 7% is larger than 
the median internal uncertainty (~ 2%), suggesting that 
either the photometric precision is dominated by uncer- 
tainties not been quantified in our data reduction or that 
many of sources contain weak excesses. We conserva- 
tively adopt a 7% uncertainty on the -R24/4.5 ratio for all 
sources. 

To identify stars that have a 24/im excess, we required 
that the residuals exceed the photospheric color by 4a 
= 0.12 dex, or a 32% excess above the photosphere. We 
adopted 4<r as the cutoff since the dispersion in the resid- 
uals can vary by ~ 0.003 dex depending on which sample 
of stars are used in the fit. By choosing a 4<r cutoff, 
weak excesses are likely significant even if the dispersion 
is underestimated by 0.003 dex. In total 54 stars have 
a value of -R24/4.5 that exceeds the stellar photosphere 
by more than 4<r and have reliable photometry. Of these 
54 stars, 35 have an 8/im and/or 16/im excess, while 19 
have a 24/iim excess only. These 19 sources are indicated 
by the filled stars in Figure [TJ The sources with excesses 
only at 24/im span the full range of spectral types and 
include 4 A/B-type stars, 1 F2 star, 3 G-type stars, 2 
K-type stars, and 9 M-type stars. 

In addition to the checks described in Section [27X1 we 
visually re- inspected the 24/im mosaics and the PRF sub- 
tracted images for the 54 stars with an infrared excess. 
Cutouts of the 24/xm images and PRF-subtracted images 



for these sources are presented in Figured) We identified 
7 stars where nebulosity or source confusion may com- 
promise the photometry and create an apparent excess. 
Three stars (HIP 77545, HIP 80024, and ScoPMS 45) 
are located in regions of extended nebulosity, although 
the observed PRF is point-like. Four other stars 
([PBB2002] USco J155729.9-225843, [PBB2002] USco 
J160532.1-193315, [PBB2002] USco J160708. 7-192733, 
and [PBB2002] USco J160900.0-190836) have a neigh- 
boring source within ~ 20". We have re-reduced the data 
for these sources using smaller sky annuli for the PRF fit, 
and examined various aperture measures. In each case, 
the infrared excess is robust to these variations in the 
data reduction. Also, in 4 of the 7 sources, the excess is 
confirmed by IRAC 8/im or IRS 16/im photometry, where 
confusion and nebulosity are less of a concern. We there- 
fore classify these 7 stars as having an 24/im excess, but 
we note that the excesses around HIP 77545, [PBB2002] 
USco J160708.7-192733, and ScoPMS 45 are detected 
only at 24/im and are otherwise unconfirmed. 

The ratio of the observed 24/im flux density to the pho- 
tospheric value was estimated from the observed photom- 
etry and the expected stellar colors. If a source does not 
have an excess in the 4.5/im or 8/im bands (see Paper I), 
the observed-to-photospheric 24/im ratio was computed 
from the observed i?24/4.5 ratio and the photospheric 
value determined from the -R24/4.5 vs. (J — -Ks)o relation 
shown in Figure [JJ Otherwise, the 24/tm photospheric 
flux density was bootstrapped from the its-band pho- 
tometry (or J-band if a K s excess is present) using em- 
pirical rel ations between spectral ty pe and J — K s color 
(Paper I, Da hm fc Carpent er 2009), spectral type and 
K — M color 6 compiled by IKenvon fc Hartmannl (|1995H , 
and the -R24/4.5 vs. (J — K S ) Q relation. Table [2] lists the 
ratio of the observed-to-photospheric 24/im flux densities 
for each star. 



3.2. 10[im excesses 

The photospheric flux density at 70/im was estimated 
from the 24/im flux density by assuming that the stellar 
emission is in the Rayleigh- Jeans limit (i.e. S v oc v 2 ). If 
the source has a known 24/im excess, the photospheric 
24/im flux density was estimated as described in Sec- 
tion [XT] Otherwise, the observed 24/im flux density was 
used. 

Figure [3] presents a histogram of the measured signal- 
to-noise ratio for the observed 70/im infrared excess 
above the stellar photosphere. The detected 70/im flux 
density from the B0. 5 star HIP 78820 is consistent with 
the photospheric emission. Excluding the two stars 
where the measured 70/im flux density is suspect (see 
Section 12. 2p . 12 stars have an apparent 70/im excess 
greater than 3a. Sources with 70/im excesses include 
4 A/B stars, 6 K stars, and 2 M stars. All 12 of these 
stars also have a 24/im excess, as well as an 8/im or 16/im 
excess identified in Paper I. 

6 We compared the empirical K — M colors versus the model K — 
[4.5] colors from the Stellar Performance Estimation Tool (STAR- 
PET) available on the Spitzer web site. Between spectral types 
of AO and G5, the maximum difference between the two colors 
is 0.03 mag. For M0 stars, the color difference is 0.15 mag. A 
comparison with the observational data indicates that the K — M 
colors provide a better match to the data. 
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4. NATURE OF THE DISK CANDIDATES 

The circumstellar dust producing the infrared excesses 
observed in Upper Sco may originate from either a "pri- 
mordial" or "debris" disk. Conceptually, primordial 
disks contain copious amounts of gas and dust that are 
remnants of the star formation process, while debris disks 
are gas-poor systems created by the collisional shattering 
of planetesimals. 

Distinguishing primordial from debris dust is difficult 
observationally at an age of 5 Myr. The inner por- 
tion of most primordial disks are optic ally thin by this 
age as traced by 3-10/im ob servations (jMamaiek et al.1 
12004 Hernandez et "all 120071 ). but the state of the outer 
disk as traced by wavelengths longer than 24/im re- 
main poorly characterized. The classic argument to sup- 
port the debris interpretation is that radiative and colli- 
sional processes can deplete the dust on timescales much 
shorte r than the stellar age ass uming the disk is gas 
poor (Backma n fc Parescd 119931) . Therefore, the dust 
must be replenished continuously, presumably from a 
collisional cascade, to explain the observed infrared ex- 
cesses. This argument assumes that gas drag is negligi- 
ble, which can circularize the orbits of the dust grains and 
reduce the frequency of destruc tive collisions. Models by 
iTakeuchi fe Artvmowic9 (|2001[ ) indicate that 10 M e of 
gas distributed radially over ~ 100 AU can prolong the 
lifetime of grains around A-type stars to a couple of mil- 
lion years, which is comparable to the age of Upper Sco. 

Unfortunately, observations of the gas component are 
not available for most of the Upper Sco sample. Several 
stars in Upper Sco have strong Ha emission lines that 
suggest gas accretion is present, although most stars do 
not appear to be accretin g detectable levels of material 
(jDahm fc Carpenteri 11)09 ]). Yet even these observations 
only trace gas in the inner disk and are not direct diag- 
nostics of the total gas surface density. Extensive surveys 
for gas in 5-100 Myr stars have been presented in other 
studies that are designed to trace the gas over a range of 
orbit a l radii (IZuckerman et al.in~995l ; INaiita fc Williams! 
120051: IPascucci et al.l 120061) . Even these observations, 
though, cannot rule out that gas drag has a significant 
influence on the dust dynamics. 

Given the lack of observations of the gas component, 
we rely on distinctions in the spectral energy distri- 
butions between stars in Upper Sco and other known 
sources to identify potential debris systems. In Paper I, 
we compared the spectral energy distributions of the Up- 
per Sco sample between 2 and 16/im with Herbig Ae/Be 
stars and classical T Tauri stars. We suggested that the 
infrared excesses around the B/A stars represent debris 
systems, while the K/M stars with 8/im and 16/im ex- 
cesses originate from optically thick disks, albeit with 
suppressed levels of mid-infrared emission relative to a 
typical classical T Tauri star in Taurus. The 24/im pho- 
tometry and IRS spectra for a subset of these sou rces 
support these conclusions (jDahm fc Carpenterll20 09). 

With the discovery of 19 additional disks in Upper Sco 
from the 24/im photometric survey, we further investi- 
gate the range of disk characteristics. Figure 0] plots the 
24/im to 4.5/im flux ratio versus the 8/im to 4.5/im flux 
ratio for stars in Upper Sco with (solid black circles) and 
without (gray circles) 24/im excesses. Qualitatively, a 
dichotomy is present in the excess characteristics. One 



population, present only for stars with spectral types 
later than K0, has both 8/im and 24/im excesses, and 
has colors similar to stars surrounded by optically th ick, 
primordial disks (Paper L lDahm fc Carpenterll2009l ). 

The second population of 24/im excess sources, present 
among all spectral types, has weak or no detectable 8/im 
excesses. Among the K0-M5 stars, the magnitude of the 
24/im excess for this second population is less than the 
24/im excess found in any of the sources with both 8/im 
and 24/im excesses. If we assume the dust is isothermal, 
a lower limit to the dust temperature obtained from the 
24/im excess and the 70/tm upper limit is ~ 53 K for 
blackbody emission. The upper limit to the fractional 
dust luminosity for this dust temperature and a stellar 
temperature of 4000 K is ~ 10~ 4 . This fractional lumi- 
nosity is orders of magnitude lower than that found from 
optically thick disks around classical T Tauri stars. 

With the large number of Spitzer surveys of 
star forming regions, the primordial disk classifi- 
cation has been refined to include "transitional" 
(iStrqm et al.l [19891 ) , "pre-transitiona l " fE spaillat ~t al.l 
|2007j), "evolved" (Hernandez et al.l I2007D . "homolo- 
gously depleted" (jCurrie et al.1 l2009a| ). and "anemic" 
(jLada et al.ll2006| ) designations. Each of these disk types 
is characterized by reduced levels of emission at wave- 
lengths < 10/iin relative to classical T Tauri stars. The 
reduced infrared emission can be caused by a number 
of processes, including dynamical clearing of dust by a 
companion, grain growth, and disk settling. 

To compare the characteristics of the Upper Sco sam- 
ple w ith anemic disks, Figure | H shows the infrared colors 
from (jLada et al.ll2006l see also lCurrie fc Kenvonll2009bD 
for c onfirmed members of IC 348 (see iLuhman et al.l 
2003). IC 348 was chosen for this comparison since the 
cluster contains a similar range of spectral types as Upper 
Sco, and was used to define anemic disks. One limita- 
tion, however, is that the IC 348 observations detected 
the photosphere at 24/im for A-type stars only, and are 
not sensitive to the weak 24/im excesses observed around 
the later type stars in Upper Sco. The vertical dashed 
lines show the range of slopes in the IRAC bands used 
to define anemic disks, while an 8/im to 4.5/im flux ratio 
great er than 0.64 would be classified as an optically thick 
disk (jLada et al.ll2006[) . None of the B/A/F/G stars in 
Upper Sco with an infrared excess would be classified 
with anemic or optically thick disk based on these crite- 
ria. The circled crosses in Figure [4] show the actual st ars 
in IC 348 that were identified bv LLada efall (|2006| ) to 
have anemic disks based on the fitted slope to all 4 IRAC 
bands. For B7-A9 (top panel in Figured]) and F/G spec- 
tral types (middle panel) , the anemic and optically thick 
disks in IC 348 have larger 8/im excesses than any of the 
sources in Upper Sco. The 24/im excesses also tend to 
be larger in IC 348, although the distribution overlaps 
significantly with Upper Sco. 

The Upper Sco sources with excesses detectable at 
24/im only also do not fit the operating definition of tran- 
sitional disks. Transitional disks are noted by a lack of 
continuum excesses at wavelengths < 10/im, but retain 
excesses at longer wavelengths comparable in strength 
to classical T Tauri stars. As shown in Figure H] the 
K/M stars in Upper Sco with excesses at 24/im only have 
systematically lower 24/im excesses than the stars with 
optically thick disks. This does not necessarily indicate 
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that the dust is debris, as a spectrum of inner hole sizes 
may exist, and the weak 24/im excesses may simply indi- 
cate a larger inner hole is present in the disks around the 
Upper Sco sample. While none of these systems were 
detected at 70/im, most primordial disks were not de- 
tected at 70/im either. Therefore, the 70/im limits for 
the weaker 24/im excesses are insufficient to rule out a 
cold, optically thick outer disk. 

In summary, the above discussion suggests that Upper 
Sco contains two populations of excess sources. The late 
type stars with strong 8/im and 24/im excesses plausibly 
have primordial disks based on the presence of accre- 
tion signatures in a few stars, and a similarity in the 
infrared colors of yo ung stars surrounded by optically 
thick disks (Paper I. IDahm fc Carpenter! I2009D . Many 
of the se sources have in frared colors similar to "anemic" 
disks (jLada et al.l l2006). The origin of the dust produc- 
ing the weak 24/im excesses remains uncertain, but these 
disks are likely not a simple extension of the transitional 
or anemic disk designations. Without additional obser- 
vations, we can only assume that these disks are gas-poor 
and the dust originates from a debris disk, but we cannot 
rule out that the dust is remnant primordial material. 

In Table [2l we classify each of the disks candidates 
as primordial, debris, or Be. Primordial disks were as- 
signed to K/M type stars that have 8/im and 16/xm ex- 
cesses since their infrared excess ch aracteristics are sim- 
ilar t o optically thick disks (see also IDahm fc Carpenterl 
[2001 Paper I). Two stars (HIP 77859 and HIP 78207) 
have prominent hydrogen emissions in opti cal spectra 
and have been identified a s classical Be stars (jCramptonl 
119681 : [Jaschek et aL[|1964D ; the observed infrared excess 
may originate from free-free e mission in an ionized gas 
disk dPorter fc Riviniusl l200l . The remaining sources 
with infrared excesses were classified as debris disks given 
the relatively weak excesses at 8, 16, or 24/im. As noted 
in Paper I, optically thick primordial disks are present 
in 19% of the K/M stars in the Upper Sco sample, and 
in none of the B, A, F, and G stars. Debris disks are 
detected around all spectral types, with percentages of 
9% for K/M stars, 13% of F/G stars, and 28% of B/A 
stars. 

5. PROPERTIES OF THE DEBRIS DISKS 

A direct comparison of the debris fractions between 
early and late type stars does not necessarily inform how 
the frequency of disks varies with spectral type. Since 
observations at a given wavelength probe dust at smaller 
orbital radii around cooler stars for a fixed grain size, the 
low 24/im debris fraction around K/M stars relative to 
A-stars may simply reflect that the debris is located at 
large orbital radii and is too cool to detect in the 24/im 
band. 

Nominally, the orbital radius of the debris can be in- 
ferred from the dust temperature. This is possible only 
for the debris disks around the B/A stars where the de- 
bris emission is detected at two or more wavelengths. For 
the F/G/K/M stars, the debris emission is detected only 
at 24/xm, and a broad range of temperatures and radii 
can fit the single photometric point. We instead pose 
the question: can the debris properties inferred around 
the B/A stars also explain the debris emission observed 
around the later spectral types? 



The emitted radiation from a debris disk varies with 
spectral type due to differences in stellar heating and ra - 
diative "blowout" of small grains (jPlavchan et al.l l2009) . 
To compare the debris properties between spectral types 
in a self-consistent manner, we assume the debris is lo- 
cated in a narrow ring at an orbital radius R. The par- 
ticles follow a power-law size distribution, n(a), between 
grain radii a m j n and a max such that 

n(a) = Nj—) a , (1) 

^ ^min ' 

where N Q is the normalization constant that effectively 
determines the total grain surface area. The power- 
law exponent is fixed at a = —3.5 as appropriate for 
an infinite collisional cascade without binary conditions 
on th e minimum and maximum particle size (Dohnanvi 
1969). In practice, the maximum grain size was fixed at 
a-max — 1000/im, and the minimum grain size (a m i n ) was 
set to either the blowout size or 0.05/iin, which ever was 
larger. Excluding smaller grains does not impact the re- 
sults significantly since such small grains are inefficient 
emitters at 24/im. 

The grain sizes that will remain gravitationally bound 
to a star were assessed by balancing the stellar gravita- 
tional force (F s ) against repulsive forces from radiation 
(Fp T ) and stellar winds (F sw ). The ratio of these forces 
{0) for a particle of radius a is given by 

at \ -^P r Fsw 
Pi ) = - J! —B 

_ 3-L«(Q ra d(a) + Qsw(a) M svr v sw c/ £*) 

WirGM^pc 1 ' 

where is the stellar luminosity, M sw is the stellar 
mass loss rate, v sw is the stellar wind velocity, Q r ad(a) 
is the particle cross section to radiation in units of the 
geometric cross section, Q sw (a) is the analogous cross 
section for stellar winds, and p is the grain density 
(IBurns et all Html IGustafsonl fl99i IStrubbe fc Chiang 
2006). Assuming the disk is optically thin, particles 
for which (3 < 0.5 remain gravitationally bound to the 
star, while particles with > 0.5 are ejected. The stel- 
lar wind velocity was set to the escape velocity (v sw = 
y/2GM*/ R*). The stellar mass loss rate is uncertain 
by orders of magnitude for these young ages across all 
spectral types. The force from stellar winds is negligi- 
ble for the grain radii considered here unless the mass 
loss rate is > 1000 times the solar value ([Plavchan et alJ 
[20051: IStrubbe fc Chiangl l200l . We assume M sw = 
2 x 10~ 14 Affn yr~\ w hich corresponds to the solar mass 
loss rate ()Woodl 120041 ). and that Q sw (a) = 1. The dust 
emission was com puted using optical constants for astro- 
nomical silicates ()Weingartner fc Draind 120011 ) and as- 
suming p = 2.7 g cm -3 . Stellar photospheres were ap- 
proximated as blackbodies. 

Given the debris disk model, we can normalize the two 
free model parameters (R and N ) using the nine B/A 
stars in Upper Sco that have both 16/iin and 24/im ex- 
cesses. The ratio of the 16/im to 24/im excess flux den- 
sities determines the orbital radius of the dust, R. The 
ratio varies between 2.8 and 10.0 with a median value 
of 3.8, which leads to an orbital radius between 9 AU 
and 40 AU with a median of 15 AU. Given the orbital 



Debris Disks in Upper Sco 



7 



radius, the total surface area of dust contained in parti- 
cles gravitationally bound to the star {[3 < 0.5) is set by 
the magnitude of the 24/im emission given the inferred 
orbital radius, which determines N a in Equation [TJ The 
24/im excess for all 13 B/A stars with debris disks varies 
between 0.5 and 16 times the photosphere, with a median 
of~2. 

To compute the excess emission around stars of other 
spectral types, stellar properties (mass, luminosity, and 
temperature) were o btained from the 5 Myr pre-main- 
sequence isochrone of lSiess et al.l (|2000t ) for solar metal- 
licity. The orbital radius (R) and the size distribution 
(N ) of the dust were fixed to that inferred around the 
B/A stars. However, for a given L* and M*, only grains 
with (3 < 0.5 were included in the calculations of the de- 
bris emission. Therefore, even though N a is fixed, the 
total dust surface area is larger for later spectral types 
since the radiation blowout size is smaller. The ratio 
of the expected dust emission excess to the photosphere 
was computed in this manner for stars between masses 
of 0.1 M© and 7M . 

Figure [5] presents the results of the calculations for 
16/im and 24/im excesses. The shape of the model curves 
reflects variations in stellar heating and the minimum 
grain size present in the disk. The dust surface area de- 
creases toward early spectral types and tends to reduce 
the dust emission, which is compensated by warmer dust 
temperatures. The 24/im excess relative to the photo- 
sphere (hereafter referred to as the relative 24/im excess) 
peaks at a spectral type of ~B7 (J — K s rj —0.06), and 
declines toward more luminous stars as more small grains 
are ejected by radiation pressure. Toward later spectral 
types, the relative 24/im excess at first declines despite 
the increased surface area since the dust temperature 
decreases with lower stellar heating. However, at a color 
of J — K s w 0.2 (~F0 star), the blowout particle ra- 
dius is ~ 5/im, and the dust temperatures are warmer 
than expected based on blackbody radiation. Combined 
with the increased surface area of these smaller grains, 
the normalized 24/im excess ratio plateaus toward cooler 
stars until a stellar color of J — K s w 0.8 (~M0 star). 
At this point the blowout size is less than 0.1/zm and 
these small grains emit inefficiently at 24/im. There- 
fore the additional surface area from small grains does 
not appreciably increase the amount of infrared excess. 
The 16/im curves exhibit similar behavior, but fall more 
steeply toward later spectral types since the dust emis- 
sion is further into the Wien regime and more sensitive 
to the dust temperature. 

Figure [5] shows that if the only variation in the de- 
bris properties is the blowout size and stellar heating, 
the 24/um excess relative to the photosphere will be ~ 2 
times higher around ~ AO stars compared to F0-M0 stars 
at an age of 5 Myr. The relative 24/im excess decreases 
rapidly toward stars earlier than B7 and later than M0. 
Qualitatively, the range of relative 24/im excesses around 
the F/G/K/M-type stars are consistent with the extrap- 
olation of the excesses around the B/A stars based on 
this model. At 16/im, the relative photometric excess 
falls sharply with later spectral types and 16/im excesses 
would be undetectable toward the later type stars. This 
model is consistent with the lack of 16/im excesses from 
debris in F/G/K/M stars (Paper I). We conclude that 
is not necessary to invoke systematic differences in the 



radial distribution of dust with spectral type to explain 
the range of observed 16/im and 24/im excesses. 

6. EVOLUTION OF DEBRIS EMISSION 

Variations in the debris luminosity with stellar age pro- 
vide a means to investigate the evolution of planetes- 
imal belt s and potentially the fo rmation of planetary 
system s. iDominik fc Decinl ((20031 see also iWvatt et al.l 
l2007b| ) suggested that the evolution of debris emission 
can be explained by quasi-steady state collisional equi- 
librium of planetesimals distributed in a narrow ring. 
One limitation of this model is that it does not pre- 
scribe a physical mechanism to initiate the coll i sional 
cascade. In a series of papers. iKenvon fc Bromlevl ( 2008. 
and references therein) presented a model that follows 
both the collisional growth of planets and the onset of 
debris production. They show that the debris disk ini- 
tially has a low luminosity when planets are still in the 
formative stages. When the collisional growth produces 
~ 2000 km size bodies, the resulting gravitational inter- 
actions excite planetesimal collisions that leads to dra- 
matically increased debris luminosity. The timescale to 
form a 2000 km body increases with the orbital radius, 
thereby producing a "wave" of debris production in time 
that propagates from the inner disk outwards. 

A prediction of the Kenyon & Bromley model is that 
the debris luminosity initially increases with time, and 
then declines once the wave of debris production has 
propagated outward through the disk on a dynamical 
timescale. For a debris disk that extends between 30 
and 150 AU, the debris emission will peak at an age 
of ~ 5-30 Myr, with faster timescales for higher mass 
stars for a fixed disk mass. A decline in debris lu- 
minosity for ages older than ~ 20 Myr has been ob- 
served in a numb er of studies and seems firmly estab- 
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further suggest that the 24/im excess around A/F stars 
increases from an age of ~ 5 Myr and peaks around 10- 
15 Myr. If verified, this initial increase in 24/im debris 
luminosity may signify the onset of planetesimal stirring 
in debris disks. 



The Upper Sco observations presented here provide 
an important data point to re-evaluate the trends at 
young ages since few surveys for debris disks are avail- 
able for younger stars than 10 Myr. As described 
in the Appendix, we have compiled the results from 
a number of Spitzer 24/im surveys for comparison to 
Upper Sco. Simil a r compilation s have appeared in 
iGasoar et all (l2009h. iRebull et all (|2008h . iCurrie et all 
(2008al). and iSiegler et all (|2007h among others. The 
compiled surveys include clusters and associations be- 
tween ages of 2 and 757 Myr, and main-sequence field 
stars between ages of 3 and 5000 Myr. The young as- 
sociations and clusters that have been surveyed previ- 
ously are eith er sparsely populated with only ~ 20 mem- 
bers ( n Cha. iGautier et all l2008t TW Hydra moving 
group, ILow et al.l 12*005) or are at large distances such 
that Spitzer is sensitive to the stella r photosphere for 
only A -type stars or earlier ( IC 348, lLada et aL 2006; 
a Ori. iHernandez et all 120071 : Ori on OBla and OBlb , 
Hernandez et al.ll2006t NGC 2262. lCurrie et al.ll2008bl) . 
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The spectral type bins were selected to be B7-A9, F0- 
F9, and G0-K5, which correspond to the spectral type 
ranges where numerous Spitzer studies exist (see Ap- 
pendix). Stars earlier than B7 were omitted since they 
are fewer in number, and as shown in Section [5J the de- 
bris luminosity may be depressed by radiative blowout 
of the smallest grains. Stars later than K5 are omitted 
since most Spitzer surveys of young clusters and associa- 
tions were insensitive to the stellar photosphere for K/M 
stars. The relationship between the empirical variable of 
spectral type and the physical variables stellar mass and 
luminosity, which dictate debris disk properties, will be 
discussed in Section 

As discussed in Section [4j we are faced with the am- 
biguity of distinguishing if the "anemic" , "evolved" , and 
"transitional" disks are debris or primordial systems. We 
proceed by an entirely observational definition. The disk 
around HR 4796A is commonly assumed to be debris in 
nature. The ratio of the observed 24/mi emission to the 
photosphere for this disk is 97, which is one of the largest 
relative 24/xm excesses known. We therefore assume that 
any B/A/F star with an "anemic" , "evolved" , or "transi- 
tional" classification and has a observed-to-photospheric 
24/im flux ratio of < 100 is a debris disk, while more lu- 
minous systems are primordial disks. For G and K-type 
stars, we use a limit of 5 based on the division shown in 
FigureU] These divisions are arbitrary and not physically 
motivated. Changing these boundaries by even a factor 
of two has no substantial impact on our conclusions. 

6.1. 24^m excess fraction 

We first examine how the fraction of stars with 24/im 
infrared excesses varies with age. We consider two differ- 
ent thresholds to identify infrared excess: 32% to incor- 
porate the largest number of studies, and 15% to increase 
sensitivity to fainter disks. The adopted thresholds for 
each survey are listed in the Appendix. Any surveys 
with sensitivity limits larger than these thresholds were 
omitted from the appropriate plot. The excess fraction 
is defined as the ratio of the number of debris disks to 
the total number of stars (i.e. optically thick disks, de- 
bris disks, Be stars, and non-excess stars). Stars were 
omitted from the statistics if the original study noted 
that the photometry was contaminated (e.g. nebulos- 
ity, nearby bright star) such that is was not possible to 
determine if the star has an excess or not. 

Figures [6] presents the excess fraction for the three 
spectral type groupings. For G0-K5 stars, the 24/xm 
excess fraction in Upper Sco is 11% (6/54) at the 32% 
detection threshold, which is indistinguishable from the 
excess fraction (6/56) observed for 10-30 Myr G0-K5 
stars. This result is consistent with the FEPS program, 
which found that the excess fraction fo r solar-type stars is 
roughly constant for a ges < 300 Myr (jMever et alJl200ll 
iCarpenter et al.ll2009ft . 

Other studies have found that the debris disk fraction 
for solar-type stars ca n be as high as 40-45% between 
ages of 3 and 30 Myr (jSiegler et al.ll2007t iGaspar et all 
2009). The difference can be attributed to semantics, 
as these studies have defined "solar-type" stars as FGK- 
stars. Whether or not F-stars should be included in the 
definition of solar-type stars is arbitrary. We treat them 
separately since there is evidence that the excess fraction 
of F-type stars is higher than G0-K5 stars. In the 10- 



30 Myr bin, the excess fraction for F-type stars is 15/37 
(41%) at the 32% detection threshold, compared to 6/56 
(11%) for G0-K5 star in the same age bin. The proba- 
bility that these two observed excess fractions are drawn 
from the same parent population is p — 0.002 according 
to the two-tailed Fisher exact test. 

For B7-A9 and F0-F9 spectral types, the 24/im excess 
fraction peaks at ~ 50% for ages of ~ 10-30 Myr and 
a 32% detection threshold. The possible peak in th e 
excess fraction result was noted bv lCurrie et al.1 <|2008bf ). 
and the declin e in the exce s s frac tion toward older ages 
was found by iRieke et alJ (l2005f). No such incr ease is 
apparent for the G/K stars IGaspar ~et~ai1 (|2009h found 
that the excess fraction for F-type stars was greater for 
ages than 3-10 Myr than for 10-30 Myr, which can be 
attributed to the high excess fr action for F-type stars i n 
Orion OBlb and Orion OBla , (jHernandez et al.ll2006h . 
Since iHernandez et al] (|2006f ) estimate that the Orion 
24/im observations are complete for stars earlier than F0 
in Orion OBla and A8 in Orion OBlb, the excess fraction 
for F-stars in Orion may be biased toward higher values 
and are not included in this analysis. 

To evaluate if the excess fraction in the <10 Myr and 
10-30 Myr age bins could have been drawn from the same 
parent population for earlier spectral types, we compared 
the ratios using the two-tailed Fisher's exact test. For 
a 24/im excess threshold of 32%, which has the largest 
sample of stars, the probability that these two age bins 
are drawn from the parent population is 17% for both B7- 
A9 and F0-F9 spectral types. If we group B7-F9 spectral 
types together, the corresponding probability is 9%. The 
significance of the enhancement in the excess fraction for 
10-30 Myr stars relative to <10 Myr stars then is ~ 1.7cr 
for stars earlier than F9. Therefore, we conclude that no 
strong evidence exists for a change in the excess fraction 
from 2-10 Myr to 10-30 Myr among early-type stars. 

6.2. 2\\im excess 

Evolution in the debris disk properties may manifest 
itself in the excess luminosity as well as the excess frac- 
tion. Figure [7] shows the ratio of the 24/xm excess to 
the photospheric flux as a function of stellar age. Only 
stars brighter than the photospheric detection limit for 
the respective samples are shown (see Appendix). Dif- 
ferent symbols are shown for primordial (asterisks) and 
debris (solid circles) disks; 3cr upper limits are shown for 
stars without a detected excess (open triangles). 

For B/A spectral types, the uppe r envelope o f 24/xm 
excess emission declines with time. IRieke et al.l (2005) 
suggested that for ages older than 5 Myr, the upper en- 
velope varies with stellar age as t , which is shown as 
the dashed line in the top panel of Figure [7] and extrap- 
olated to younger ages. The relative 24/im excesses ob- 
served around B/A stars in Upper Sco are consistent 
qualitatively with this upper envelope. For G0-K5 spec- 
tral types, the peak relative 24/im excess is roughly con- 
stant for ages less than 300 Myr, with a possible decline 
toward older ages. For F-type stars, a sharp peak in the 
24/im excess i n present at an age of ~ 16 Myr as previ- 
ously n oted bvfHcrnand ez et al.1 (|2006l ) and lCurrie et al.1 
( 2008a). This peak is due to lumino us debris disks in t he 
Lower Centaurus Crux association (|Chen et al.l l2005). 

To evaluate the trends quantitatively, we used the 
generalized Kendall's Tau test as implemented in the 
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ASURV Rev 1.2 package (jLavallev et al.l fl993) . This 
test incorporates both the detections and upper limits 
to make efficient use of the data. We first consider if 
evolution is present over all ages in the sample. The 
probability that a correlation of the 24/im excess versus 
age is not present for B7-A9, F0-F9, and G0-K5 stars 
is 3 x 10~ 10 , 2 x 1(T 8 , and 1(T 5 , respectively. These 
results confirm the trend apparent from visual inspec- 
tion of Figures [7] that the magnitude of the 24/zm excess 
decreases with stellar age over all spectral types consid- 
ered here. Th is conclusion has been reached previousl y 
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We now consider if evolution in the 24/im excess is 
present at young ages. We specifically consider the 
age range between 5 and 17 Myr, which is well pop- 
ulated with clusters and associations, and significant 
evolution is anticipated based on theoretical models 
for the collisional evolu tion of a Kuiper Belt analog 
(jKenvon fc Bromlevi[200l) . For B7-A9, F0-F9, and G0- 
K5 spectral types, the probability that a correlation of 
the 24/itm excess versus age is not present between 5 
and 17 Myr is 0.14, 0.01, and 0.56, respectively. Thus 
when subdivided by spectral type, the only suggestion 
of a trend in the 24/im excess with age between 5 and 
17 Myr is for F-type stars. From inspection of Figure [3 
the trend is such that the amount of 24/zm excess in- 
c reases over this age ra nge . 

iCurrie et al.1 (|2008aD reported a significant rise in the 
magnitude of the debris emission between 5 and 17 Myr 
based on debris disk observations around B/A/F stars 
in Orion OBlb, Orion OBla, LCC, and UCL. They used 
the Wilcoxon rank-sum test to determine the probabil- 
ity (p) that the mean excess is consistent between two 
samples. The stron gest statistical trend identified by 
ICurrie et all C2008a) was the increase in the mean 24/im 
excess around A/F stars in the 5 Myr Orion OBlb asso- 
ciation compared to the 8.5 Myr Orion OBla association 
(p = 0.002), and Orion OBlb compared to a combined 
16-17 Myr old LCC/UCL sample ( p = 0.05). One differ- 
ence in the analysis conducted by ICurrie et al.1 (|2008a) 
and this study is that we separated the F-type stars from 
the B/A stars, and we did not include the F-type stars in 
Orion OBla and Orion OBlb. While the separation by 
spectral type is arbitrary from a physical point of view, 
the observations are complete for B/A spectral types over 
a broader range of ages. The MIPS 24/im observations of 
Orion OBla and Orion OBlb in partic ular are complete 
to th e photosphere for F0 and A8 stars (| Hernandez et all 
2006), respectively, and the 24/zm detections of F-type 
stars in these associations may be biased toward stars 
with disks. 

We stress that we have not included in this analysis the 
17 debris disks identified in the MIPS 24/im survey of the 
13 My r old, double cluster h & x Persei (jCurrie et al.1 
I2008a[) . These 17 stars have A/F spectral types with 
relative 24/im excesses that range between 8 and 180, 
with a median value of 19. These excesses are extreme 
compared to the sources shown in Figures [7J These lu- 
minous debris disks could conceivably indicate a peak in 
the debris production in the 10-20 Myr age range. How- 
ever, the parent sample for the h & x Persei is about 



30 times larger than the Upper Sco sample. As noted by 
ICurrie e~ al. (2008a|), the extreme excesses in h & x Per- 
sei may simply reflect better sampling of the same parent 
luminosity function. Deeper observations and tabulation 
of the measured flux densities for all cluster members will 
determine if the extreme excess sources in h & x Persei 
represent a peak in the debris production. 

6.3. Interpretation 

Before drawing further conclusions from the results 
presented in Section 16.21 wc consider how the stellar 
properties may influence the interpretation of debris disk 
evolution. The data presented in Figures [HI and [7] were 
grouped by the observed spectral type. For main se- 
quence stars, the spectral type bins of B7-A9, F0-F9, and 
G0-K5 correspond approximately to stellar masses of 1.9, 
1.4, and 0.8 M Q . The corresponding Kelvin-Helmholtz 
contraction times are ~3, 10, and 30 Myr. Since the 
ages of the stellar samples are as young as 2 Myr, the 
younger stars are in the pre-main-sequence phase of evo- 
lution and the stellar mass and luminosity will vary with 
age for each of the spectral type bins. This is particu- 
larly true for young (< 10 Myr) stars more massive than 
~ 1 M , which are evolving toward the main-sequence 
along the Henyey tracks. For example, the spectral 
type for a 2M Q star will be K2, G8, and A2 at an 
age of 2, 5 and 10 Myr, respectively according to the 
ISiess et al.1 (12 000) pre-main-sequenc e evolutionary mod- 
els and lKenvon fc H artmann (1995) temperature scale. 

Since the blowout size is proportional to L*/M* (see 
Equation^, the minimum grain size in the debris disk 
will also vary with age. For a fixed spectral type the 
I/*/M* ratio, and hence the blowout particle size, de- 
creases by a factor of ~ 5 between 1 Myr and 10 Myr 
for A0-M0 stars. The total grain surface area is pro- 
portional to a~? n 5 for a a -3 5 particle size distribution, 
and the debris luminosity will thus increase in time if the 
planetesimal belt is in otherwise steady state. 

The expected trends were evaluated more quantita- 
tively using the debris disk model described in Section [5j 
For a given age and spectral type (i.e. stellar tempera- 
ture), th e stellar mass and luminosity were determined 
from the ISiess et alj (|2000D evolutionary models for solar 
metallicity and no convective overshoot. The 24/im de- 
bris emission was then computed for a planetesimal belt 
at an orbital radius of 15 AU, with a minimum grain size 
corresponding to the blowout size. The results are insen- 
sitive to the assumed orbital radius since the variations 
in the debris luminosity with time are caused primarily 
by changes in the grain surface area. We assume that 
the minimum grain size in the debris disk varies instan- 
taneously as stars evolve since the collisional lifetime of 
dust within an orbital rad ius of 30 AU is less than 1 Myr 
for observed debris disks (jBackman fc Parescdl 19931 ). In 
this manner, the variation in the 24/im excess relative to 
the photosphere was computed considering only stellar 
evolution. 

Figure [5] shows the results of the model calculations. 
For AO through K0 stars, the relative 24/zm excess in- 
creases by a factor of ~ 5 from an age of 1 Myr to a 
peak at 10-50 Myr. The peak relative excess is reached 
sooner around earlier spectral types since the stars reach 
the main-sequence on shorter time scales. M-type stars 
obtain a peak relative excess in only a few million years 
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when the grain blowout size is less than 0.05/im. At 
this point, any grains removed by radiation pressure are 
inefficient emitters at 24/im, and the debris luminosity 
declines with the stellar luminosity. 7 If the model as- 
sumptions are correct, the sensitivity functions presented 
in Figure [8] will be superimposed on the evolution of the 
debris luminosity and must be considered when attribut- 
ing any increase in the debris luminosity with age to 
planet formation. 

For comparison, the gray solid line in Figure [8] shows 
the variation in the 24/im excess relative to the pho- 
tosphere for a 2 M Q star (~ A3 spectral type) with a 
planetesimal belt betwee n 30 and 150 AU from the 
iKenvon fc Bromlevl (|2008f ) models. Since the stellar lu- 
minosity is held constant in these models, the temporal 
variation in the infrared luminosity results from changes 
in the debris production. The debris evolution in these 
models depends on the stellar mass and the planetesimal 
belt properties (mass, orbital radius), but the general 
prediction is that the fractional debris luminosity rises at 
early ages, and declines at later times. The results show 
that for a 30-150 AU planetesimal belt, the rise in the 
debris luminosity and the age where debris production 
is maximized resembles the evolutionary curve expected 
based on stellar evolution alone. Clear differences are ob- 
served at older ages wher e the debris luminosity in the 
IKenvon &: Bromlevl (|2008f ) models declines more rapidly. 

We now evaluate if the magnitude of the relative 24/im 
excess varies with stellar age for ages < 20 Myr, and if 
any increase can be attributed to stellar evolution. For 
both the B7-A9 and G0-K5 spectral type ranges, any 
variations in the 24/xm excess with age between ages of 5 
and 17 Myr have a significance of less than 2a (see Sec- 
tion !6.2l and Figure[7]). Thus independent of the selection 
function presented in Figure El no compelling evidence 
exists for a variation in the mean relative 24/im excess 
for B7-A9 and G0-K5 stars in that age range. 

The F-type stars show the strongest evidence ( 2.6a ) 
for a rise in the relative 24/xm excess (see Section [6~2]) . 
but these results need to be interpreted in view of the ex- 
pected increase in the debris emission shown in Figure [8l 
For the F-type stars between ages of 5 and 17 Myr, we 
scaled the magnitude of the 24/im excess by the inverse 
of the model calculations computed for a F5 star, which 
corrects for the reduced excess emission expected to oc- 
cur with stellar evolution. The probability of a trend in 
the observed excess emission with age (see Figure [7]) was 
then recomputed using the Kendall Tau test. We find 
the probability that a correlation is not present increases 
from p = 0.01 (2.6a) to p = 0.02 (2.3ct). 

In summary, we conclude that any evidence for an in- 
crease in the 24/im debris emission with age between 5 
and 17 Myr rests primarily with F-type stars, and the sig- 
nificance is < 2.6tr. Any variations in the relative 24/im 
excess over this age range for B/A and G/K are signifi- 
cant at less than 2a confidence. Given the low statistical 
significance and that the increasing relative excess is ob- 
served in a narrow spectral type range, we hesitate to 

7 MIPS observations have yielded a smal l number of debri s 
disks around young, low mass T Tauri stars {Cieza ct al. 120071) . 
The low excess needs to be interpreted in context of these model 
calculations, which predict the 24/im excess will be 5 times weaker 
than a 10 Myr star if the planetesimal belt is otherwise the same. 



ascribe the apparent rise in the mean relative 24/im ex- 
cess among F-stars to the onset of planetesimal stirring. 

6.4. Debris disks around M-stars 

Evolution of debris emission around M-stars is more 
difficult to quantify since most MIPS surveys are in- 
sensitive to the stellar photosphere for these spectral 
types. Nonetheless, for completeness, we present in 
Figure [9] the 24/im excess relative to photo sphere as a 
funct ion of age f or M0-M5 stars in IC 3 48 (jLada et al.l 
120061) . a Ori ([Hernandez et "all 120071) . Upper Sco, 
■ n Cha dGautier et alj 120081 ), TW Hydra moving group 
(ILow et al.1120051). NGC 254 7 (jForbrich et al]l2008h . and 
IC 2391 (|Siegler et al.ll2007h . M-stars with upper limits 
to the 24/im excess are not shown for clarity, and because 
most studies have not reported upper limits if a star is 
not detected. 

Figure [9] again illustrates the dichotomy at young ages 
between classified primordial disks (asterisks) and debris 
disks (filled circles). The 24/im excess around Upper Sco 
M-stars is comparable to many of the stars in NG C 2547 
at an age of -38 Myr (jNavlor fc Jeffries! I2006D . How- 
ever, 5 M-stars in NGC 2547 have a 24/im excess of — 10 
times the photosphere : three of these 5 s tars have been 
observed with IRAC (jYoung et al.l l2004)l . and one has 
an 8/xm excess. By contrast, all of the Upper Sco sources 
classified with a primordial disk have an 8/im excess. It 
would be unusual for an optically thick disk to persist 
for — 38 Myr, and the debris disk is a reasonable in- 
terpretation for the NGC 2547 sources. These results 
suggest that distinguishing primordial disks from debris 
disks for M-type stars merits a more detailed examina- 
tion (cf. lErcolano et al.l[2009h . 

7. SUMMARY 

We have presented the results of a Spitzer 24/im and 
70/im photometric survey of 205 members of the Up- 
per Scorpius OB association. These data were combined 
with MIPS photometry for 15 Upper Sco sources ob- 
served by the FEPS Legacy program to provide a census 
of circumstellar disks around stars with spectral ranges 
ranking from BO to M 5 at an association age of 5 Myr 
(|Preibisch et al.ll2002D . By analyzing the 24/im photom- 
etry with 4.5/im and 8/im measurements presented in 
Paper I, we identify 54 stars that have observed 24/im 
emission that exceeds the expected stellar photospheric 
emission by 32% or more. Similarly, 12 stars were iden- 
tified with a > 3a photometric excess at 70/im; all 12 of 
these stars have a detectable 24/im excess. 

The nature of the excess sources were established based 
on the color and magnitude of the excess emission. We 
find a dichotomy of excess characteristics. One popu- 
lation, found only around the K and M-type stars, has 
strong excess emission at both 8/im and 24/im that is 
comparable to known optically thick circumstellar disks. 
A second population, found around all spectral types, has 
weak 8/im and 24/im excesses compared to op tically thick 
and "anemic" disks in the IC 348 cluster (jLada et alJ 
2006). We suggest that these weak excesses originate 
from debris disks formed from the collisional shattering 
of planetesimals, although we cannot exclude the pos- 
sibility that these systems are the remnants of optically 
thick disks. Of the 54 excess stars, we attribute the 24/im 



Debris Disks in Upper Sco 



11 



excess emission to 2 Be stars, 24 primordial disks, and 
28 debris disks. The debris disks include 11 K/M stars. 

The debris disks were analyzed to investigate whether 
the orbital radius of the presumed planetesimal belts 
vary systematically with spectral type. We modeled 
the emission with a pow e r-law distribution of grain sizes 
following the iDohnanvil (|1969f ) equilibrium distribution 
(n(a) oc a -3 5 ). The orbital radius (R) and grain size 
distribution (JV ; see Equation [T]) were set to reproduce 
the observed excess emission around nine B/A stars with 
infrared excesses at both 16//m and 24/im. The expected 
emission around stars of other spectral types were then 
estimated after allowing for variations in the stellar heat- 
ing and radiation blowout of small grains. We find that 
the magnitude of the 24/im emission observed around 
later stars can be produced by this model, indicating it 
is not necessary to invoke a different radial distribution. 
This model predicts a steep fall off in the relative 24/im 
excess for stars earlier than ~ B7 due to blowout of the 
small grains, and later than ~ MO due to reduced stellar 
heating. 

The Upper Sco results are combined with other Spitzer 
24/im surveys in the literature to reassess the evolution 
of 24/im debris emission. After consideration of both 
sample sizes and detection limits, we find a decline in the 
magnitude of the 24/xm excess relative to the photosphere 
for spectral types between B7 and K5 and ages between 
10 Myr and 5 Gyr, as has been noted in p r evious stud- 
ies (jHabing et all 120011 : iDecin et all l2003t iRieke et all 



120051: iCarpenter et alJl2009h . We also investigated if the 
24/im excess increases with stellar age for ages < 20 Myr. 
Such an increase may be indicative of the onset of plan- 
etesimal stirring |Kenyoir&2Iomle32Q08|) . We caution, 
however, that a similar increase in the excess emission 
may also result from pre-main sequence evolution; for a 
fixed spectral type, the L*/M* ratio will decrease with 
age, which leads to a higher debris luminosity as fewer 
grains are ejected by radiative forces. The observed mean 
24/im excess around F-type stars increases between ages 
of 5 an d 17 M yr as previously found by iCurrie et al.l 
(2008a], see also iHernandez et al.l l2006h . but with a sig- 
nificance of ~ 2.3-2.6cr confidence. The variations in the 
mean 24/itm excess around B/A and G/K stars over the 
same age range are significant at less than 2a confidence. 
We conclude that the observational data do not yet re- 
quire a physical mechanism to produce a peak in the 
observed debris emission. 
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APPENDIX 

COMPILATION OF Spitzer DEBRIS DISK SURVEYS 

To compare the Upper Sco results with previous observations, we compiled data from published Spitzer 24/im surveys 
for debris disks. We considered studies in which the parent sample is not biased for or against the presence of the 
debris disk. This requirement implies that the parent sample must have been selected by photospheric indicators, and 
that the Spitzer observations were sensitive to the stellar photosphere. In practice this is not strictly true for nearby 
moving group such as Pic, 77 Cha, and TW Hydra, where the moving group was initially identified based on the 
presence of an apparently isolated sources with an prominent infrared excess. We nonetheless include these groups 
since the subsequent surveys for moving group members were not biased toward circumstellar disks. 

Table [3] summarizes the surveys that were compiled and contain 5 stars or more. We include in this table the name of 
the region, the adopted age, the earliest spectral type in the survey, the latest spectral type at which the photosphere 
can be reliable detected as discussed above, the number of stars in this spectral type range, and the minimum 24/im 
excess above the stellar photosphere that could be detected. In some regions (e.g. NGC 2232, Pleiades, Praesepe), 
not all stars have spectral types available from spectroscopic observations, and the spectral type was estimated from 
optical or near-i nfrared colors . Field stars with ages between 3 Myr and 5 Gyr were compiled from the F EPS survey of 
solar-t ype stars (jMever et al.ll2006l) . th e survey of solar-ty pe fi eld stars bylBeichman et~a l. (2006) and lTrilling et al.1 
( 2008), and the survey of B/A stars by IRieke et all (|2005l ) and lSu et all (|2006l ). 

The photometry was extracted from the publications listed in Table H] No attempt was made to re-reduce the 
data under a common data reduction method or to re-assess cluster or association membership. Several studies noted 
sources where the photometry was contaminated by nebulosity or a companion sources. If these stars were excluded 
from the analysis in the original study, they were also excluded in our analysis. After compiling the source lists, we 
ran all source names through SIMBAD to identify duplicate sources appearing under different names. If photometry 
appear for a star in more than one survey, we adopted the photometry from the survey with the higher signal to noise. 

Most Spitzer studies have identified sources with excesses based on the presence of unusually red infrared colors (e.g. 
[K]-[24], [8]- [24]). Each study has a minimum detectable excess that depends on the depth of the survey, the presence 
of nebulosity in sky background, the distance to the stellar sample, and the intrinsic brightness of the star. In nearly 
all cases, we adopted the limiting thresholds adopted in the original study to identify 24/im excesses, which is normally 
taken to be three or four times the RMS scatter around the main locus of data points in a color-color diagram or a 
stellar color that is thought to represent photospheric colors (see, e.g., Figure [1]). We then identified the latest spectral 
type in which the internal uncertainty in the measured photospheric flux was less than this RMS scatter. This criteria 
removed sources where the photometric uncertainty was comparable to the excess being measured. If the threshold to 
identify an excess (= A) was expressed in magnitudes, the adopted threshold was estimated as 10 A / 25 . 
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The ages for clusters and associations were adopted from the literature based on a number of techniques, including 
pre-main-sequence evolutionary tracks, lithium depletion, and the main sequence turnoff point. The age uncertainties 
were adopted from the literature, but if no uncertainty was quoted, we assumed an age uncertainty of 0.3 dex in the 
logarithm of the age. The age uncertainties for the clusters and associations do not reflect any systematic uncertainty 
associated with the particular method. The ages for field stars were adopted from the studies presenting the Spitzer 
data. The references for the ages are provided in Table HI 

A few studies are not included in this comparison since the observations could not be analyzed in a consistent manner. 
First, we did not includes surveys where the observations were not sensitive to the photosphere for spectral types A 
or lat er, either beca use the regions are m ore distant than ~ 900 pc (Trum pler 37 and NGC 7160. ISicilia-Aguilar et all 
120061: h & x Persei. ICurrie et alJl2008at NGC 2262. ICurrie et al.ll2008b| ). or the observations were exclusively of late 
spectral types (A Orionis. iBarrado v Navascues et al.ll2007h . MIPS 24/im observations of the 7 Velorum cluster were 
not incorporated; while the obser vations were sensitive to the photosphere for stars early than ~ F8, the contamination 
from field stars is non-negligible (Hernan dez et al1l2008f ). W e also did not analyze MIPS observatio ns of weak- lined T 
Tauri stars since the samples excluded optically thick disks ([Padgett et al.l [20061 : iCieza et al.ll2007f) . 

Given the heterogeneous nature of the analyzes to identify Spitzer excesses, we occasionally had to re-analyze data 
from the original study to provide a uniform analysis. The following discussion summarizes the regions where we 
re-interpreted the observations. 

a Orionis: From the photometry presented by iHernandez et alJ (|2007f ). 97% (33 out of 34) stars brighter than 
J=10.2 mag were detected at 24/im, while only 4 of 7 stars between J=10.2 and 10.5 mag were detected. The 
completeness limit was then defined as J=10.2. For stars fainter than J=10.2, the bluest star has a dereddened color 
of J — H = 0.15. Therefore, we assumed the survey is complete for J — H <= 0.15, corresponding to spectral type 
earlier than F0. 

IC 348: Of the eight A-type stars in IC 348 studied bv lLada et al.l (|2006l ). six were detected by at 24/im, one was 
not observed at 24/im, and one has an upper limit from locally bright background. Only 57% of the FG stars and less 
than 33% of the KM stars in IC 348 have 24/im measurements, and therefore only the A-stars i n IC 3 48 were included. 
We adopted a detection threshold of 30% for a 24/xm excess based on Figure 7 in lLada et al.l (|2006[ ). 

NGC 2232: ICurrie et all ([2008bD estimated a 5cr and 10cr detection limit of 10.5 mag and 9.75 mag for their 24/im 
survey. They also find that the number of sources detected at 24/im falls rapidly for sources fainter than [24]=10.5. 
We therefore adopted [24] = 10 mag as a limiting magnitude, and find that 99% of stars brighter than J=10 mag have 
a MIPS 24/im detection brighter than this limit. For a distance of 340 pc, an age of 25 Myr, and a color excess of 
E(B — V)=0.07, J=10 corresponds roughly to a A8 star, which we adopt as the completeness limit of the survey. The 
RMS scatter in the K- [24] color fo r stars brighter than the completeness limit and not containing an infrared excess 
(as defined in ICurrie et"atll2008bD is 0.09 m ag, w hich was adopted as the photometric uncertainty. 

Praesepe: From Figure 6 in lGaspar et a l. (2009), we find that the number of stars with unphysical blue stellar colors 
increases for stars redder than r — K s =\.h. We consider only stars bluer than this limit, which corresponds to stars 
earlier than G5. 
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TABLE 1 

MIPS 24 AND 70/iM PHOTOMETERY 



Source MIPS 24^m 1 MIPS 70fj.ni 2 AORKEY Flags 





S v (mjy) 


a (mjy) 


S v (mjy) 


a (mjy) 




HD 142987 


7.89 


0.07 


5.2 


11.5 


13908224 


HD 147810 


9.67 


0.09 


-8.7 


10.2 


13906432 


HD 149598 


5.56 


0.05 


11.0 


21.6 


13909504 


HIP 76071 


10.79 


0.10 


-11.4 


8.8 


13902848 


HIP 76310 


165.40 


1.50 


373.1 


11.3 


13904896 


HIP 76633 


9.26 


0.08 


12.0 


9.0 


13907968 


HIP 77457 


8.63 


0.08 


14.7 


12.2 


13906688 


HIP 77545 


9.48 


0.09 


59.3 


58.5 


13899264 


HIP 77635 


79.15 


0.71 


-18.1 


28.3 


13899264 


HIP 77815 


10.85 


0.10 


0.0 


12.7 


13901568 


HIP 77840 


68.72 


0.62 


-58.5 


25.9 


13899264 


HIP 77858 


46.93 


0.42 


-5.1 


23.5 


13899776 


HIP 77859 


91.89 


0.83 


-18.4 


15.2 


13899776 


HIP 77900 


19.55 


0.18 


-12.5 


10.2 


13901056 


HIP 77909 


25.76 


0.23 


4.2 


15.9 


13899264 


HIP 77911 


142.60 


1.30 


253.1 


10.5 


13901568 


HIP 77960 


8.49 


0.08 


-15.6 


9.5 


13901568 


HIP 78099 


12.98 


0.12 


-21.8 


10.0 


13901568 


HIP 78104 


105.50 


0.90 


0.5 


9.7 


13898752 


HIP 78168 


34.18 


0.31 


20.4 


13.4 


13900288 


HIP 78196 


10.05 


0.09 


-1.4 


7.3 


13903616 


HIP 78207 


766.40 


7.30 


233.3 


25.4 


13899008 


HIP 78233 


6.22 


0.06 


-2.0 


8.8 


13900288 


HIP 78246 


35.33 


0.32 


15.6 


11.9 


13900032 


HIP 78265 


277.90 


3.10 


5.2 


88.2 


13897472 


HIP 78483 


8.04 


0.07 


7.1 


15.7 


13897472 


HIP 78494 


11.74 


0.11 


-13.8 


11.7 


13900032 


HIP 78530 


13.10 


0.12 


2.2 


17.9 


13902336 


HIP 78549 


12.20 


0.11 


-9.6 


12.9 


13902336 


HIP 78702 


11.90 


0.11 


2.4 


9.6 


13905408 


HIP 78809 


8.54 


0.08 


-9.3 


13.1 


13900544 


HIP 78820 


451.70 


4.10 


46.5 


13.4 


13896960 


HIP 78847 


10.42 


0.09 


5.7 


10.7 


13902336 


HIP 78877 


31.33 


0.28 


-3.0 


24.8 


13900544 


HIP 78933 


144.80 


1.30 


-4.0 


19.4 


13898240 


HIP 78956 


14.92 


0.13 


1.0 


10.8 


13905152 


HIP 78963 


8.85 


0.08 


14.7 


14.8 


13907200 


HIP 78968 


14.72 


0.13 


-2.2 


11.7 


13905152 


HIP 78996 


35.87 


0.32 


9.4 


10.7 


13900544 


HIP 79083 


17.87 


0.16 


-3.8 


9.7 


13898240 


HIP 79097 


10.20 


0.09 


-6.0 


9.9 


13906176 


HIP 79124 


12.64 


0.11 


0.7 


15.3 


13903104 


HIP 79156 


22.62 


0.20 


5.5 


9.4 


13897984 


HIP 79250 


12.02 


0.11 


-25.2 


12.0 


13905920 


HIP 79366 


9.27 


0.08 


2.6 


13.2 


13905920 


HIP 79374 


212.20 


1.90 


12.0 


21.3 


13897984 


HIP 79392 


7.28 


0.07 


2.4 


12.5 


13900800 


HIP 79404 


64.91 


0.58 


-48.1 


21.4 


13899520 


HIP 79410 


41.25 


0.37 


8.1 


11.9 


13897984 


HIP 79439 


23.86 


0.21 


32.8 


15.4 


13897984 


HIP 79462 


15.11 


0.14 


-3.8 


10.6 


13905920 


HIP 79530 


24.74 


0.22 


3.5 


14.0 


13900800 


HIP 79606 


11.48 


0.10 


-12.0 


14.3 


13904128 


HIP 79643 


9.30 


0.08 


0.5 


10.2 


13900800 


HIP 79644 


3.03 


0.04 


-16.7 


9.8 


13901312 


HIP 79733 


5.09 


0.05 


-2.0 


9.9 


13899520 


HIP 79739 


11.47 


0.10 


-17.7 


22.3 
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1 QQOO^AA 
10yUU044 


PZ99 


J160354.9-203137 




o.uy 


17 
U. 1 I 


fi 
-U.D 


12 1 


1 ,: tQ0Q7fi0 

ioyuy / uu 


PZ99 


J160357.6-203105 




Q9« no 


3 on 
o.uu 


1 87 & 
lo I .O 


19 7 

1Z. ( 


1 3Q0Q7fi0 
ioyuy i uu 


PZ99 


J160421. 7-213028 




1 (V7 ^o 

ID / .OU 


1 ^O 
1 . OU 


1 Q1 7 A 
iy 1 / .4 


9A £ 
Z4.o 


1 QQ09QQfi 

loyuzooo 


PZ99 


J160539. 1-215230 




1 Q7 
l.O f 


U.U4 


-Z. 1 


in /I 
1Z.4 


loyUZooD 


PZ99 


J160612. 5-203647 




Z.OO 


n O/i 

U.U4 


O K 

u.o 


in o 

lU.y 


1 QQ0C9/10 

iooyoZ4u 


PZ99 


J160654.4-241610 




O Oft 
Z.OO 


U.U4 


n o 

u.» 


1U.Z 


loyUU044 


PZ99 


J160703.9-191132 




O A Q 

Z.4o 


U.U4 


ft o 
-D.Z 


y.o 


1 QonQi ri/i 
loyUolU4 


PZ99 


J160831. 4-180241 




2.87 


0.04 


8.6 


14.6 


loylU (o4 


PZ99 


J160856. 7-203346 




3.35 


0.04 


-11.7 


10.6 


1 o o n o o a c\ 

13898240 


PZ99 


J160930.3-210459 




3.06 


0.04 


19.4 


13.0 


13911040 


PZ99 


J161302. 7-225744 




3.79 


0.04 


-7.3 


10.0 


13905920 


PZ99 


J161933. 9-222828 




3.63 


0.04 


-8.7 


10.7 


13903360 


RX J1548. 0-2908 




3.11 


0.04 


-4.2 


13.3 


13906688 


RX J1550.0-2312 




2.91 


0.04 


0.9 


17.7 


13911296 


RX J1550.9-2534 




5.02 


0.04 


4.8 


18.7 


13899264 


RX J1554.0-2920 




3.17 


0.04 


3.6 


9.5 


13898752 


RX J1558.1-2405A 




2.26 


0.04 


-46.3 


37.6 


13900032 


RX J1600.7-2343 




0.95 


0.04 


-9.6 


13.1 


13909248 


RX J1602 


8-2401A 




6.92 


0.06 


1.6 


18.0 


13900544 


RX J1603 


6-2245 




3.99 


0.04 


4.2 


12.5 


13902336 


SAO 183706 




10.87 


0.10 


-17.8 


25.0 


13905664 


ScoPMS 13 




3.08 


0.04 


-4.5 


14.8 


13899776 


ScoPMS 17 




43.36 


0.39 


59.0 


14.3 


13901568 


ScoPMS 23 




5.20 


0.05 


27.3 


12.8 


13902336 


ScoPMS 28 




1.49 


0.04 


1.9 


11.0 


13905408 


ScoPMS 29 




2.48 


0.04 


-2.8 


9.6 


13905408 


ScoPMS 31 




334.60 


3.00 


296.7 


11.5 


13905408 


ScoPMS 32 




0.98 


0.04 


-1.1 


11.2 


13909760 


ScoPMS 45 




5.48 


0.10 


29.0 


20.9 


13910784 



5.7 



MIPS 24/^m photometry measured using a flux calibration factor of 0.0447 MJy/sr. 
MIPS 70/im photometry measured using a flux calibration factor of 702 MJy/sr. 

Source names with the prefix [PZ99] and [PBB2002] arc from Prcibisch &z Zinncckcr (1999) and Prcibisch ct al. (2002), respectively. 
24/iUl photometry is unreliable due to bright, extended nebulosity. 
24/iin photometry is unreliable due to nearby star. 
Source is extended at 24/j,m 

70/i-m photometry is unreliable due to contaminating source or nebulosity. 
Contaminating sourcc(s) subtracted before 70/j-m photometry was measured. 
70/i-m photometry measured on AOR 13903104. 
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TABLE 2 
Normalized Flux Ratios 





SpT 


Ref 

Op 1 


Air 




70//TT1 


rr/70;/m \ 


H/AL.CBO ; 


TVrt 


HD 142361 


G3V 


2 


0.2 


1.02 


8.7 


13.9 






HD 142987 


G4 


4 


1.7 


0.99 


6.2 


13.6 






HD 146516 


GOIV 


8 


0.6 


0.99 


6.5 


18.3 






HD 147810 


Gl 


4 


1.0 


1.00 


-8.3 


9.8 






HD 149598 


GO 


4 


0.7 


1.13 


18.5 


36.2 






HIP 76071 


B9V 


2 


0.3 


0.98 


-9.9 


7.6 






HIP 76310 


A0V 


2 


0.1 


17.17 


361.0 


37.6 


Y Y 


debris 


HIP 76633 


B9V 


2 


0.3 


1.30 


12.1 


9.1 






HIP 77457 


A7IV 


1 


0.1 


0.96 


15.8 


13.3 






HIP 77545 


A2/3V 


2 


1.2 


1.79 


104.5 


103.6 


Y 


debris 


HIP 77635 


B1.5Vn 


2 


0.5 


0.99 


-2.1 


3.3 






HIP 77815 


A5V 


2 


0.7 


1.04 


0.0 


10.9 






HIP 77840 


B2.5Vn 


2 


0.4 


0.88 


-7.9 


3.6 






HIP 77858 


B5V 


2 


0.4 


0.97 


-1.0 


4.7 






HIP 77859 


B2V 


2 


0.6 


1.68 


-3.1 


2.6 


Y 


Be 


HIP 77900 


B7V 


1 


0.2 


0.98 


-6.0 


4.9 






HIP 77909 


B8III/IV 


2 


0.1 


0.99 


1.5 


5.8 






HIP 77911 


B9V 


2 


0.3 


9.64 


159.4 


17.2 


Y Y 


debris 


HIP 77960 


A4IV/V 


2 


0.6 


1.02 


-17.2 


10.6 






HIP 78099 


A0V 


2 


0.5 


1.31 


-15.7 


7.4 






HIP 78104 


B2IV/V 


1 


0.1 


0.99 


0.0 


0.9 






HIP 78168 


B3V 


2 


0.6 


0.99 


5.6 


3.7 






HIP 78196 


A0V 


1 


0.0 


0.96 


-1.3 


6.8 






HIP 78207 


B8Ia/Iab 


2 


0.1 


4.36 


12.4 


1.8 


Y Y 


Be 


HIP 78233 


F2/3IV/V 


2 


0.3 


0.95 


-2.9 


13.1 






HIP 78246 


B5V 


2 


0.2 


0.99 


4.1 


3.2 






HIP 78265 


B1V+B2V 


2 


0.2 


1.07 


0.2 


3.0 






HIP 78483 


GOV 


1 


0.4 


0.98 


8.3 


18.3 






HIP 78494 


A2mA7-F2 


2 


0.9 


0.94 


-11.0 


9.3 






HIP 78530 


B9V 


2 


0.5 


1.07 


1.5 


12.7 






HIP 78549 


B9.5V 


2 


0.4 


1.03 


-7.3 


9.9 






HIP 78702 


B9V 


2 


0.7 


1.26 


1.8 


7.6 






HIP 78809 


B9V 


2 


0.4 


1.02 


-10.2 


14.4 






HIP 78820 


B0.5V 


2 


0.7 


1.06 


0.9 


0.3 






HIP 78847 


AOV 


2 


0.5 


1.00 


5.1 


9.6 






HIP 78877 


B8V 


2 


0.1 


1.03 


-0.9 


7.4 






HIP 78933 


B1V 


2 


0.4 


0.91 


-0.3 


1.3 






HIP 78956 


B9.5V 


2 


0.6 


1.55 


1.0 


10.4 


Y 


debris 


HIP 78963 


A9V 


1 


0.4 


0.96 


15.4 


15.7 






HIP 78968 


B9V 


2 


0.5 


1.84 


-2.6 


13.6 


Y 


debris 


HIP 78996 


A9V 


2 


0.4 


3.71 


9.1 


10.4 


Y 


debris 


HIP 79083 


F3V 


2 


0.9 


1.03 


-2.0 


5.0 






HIP 79097 


F3V 


2 


0.3 


1.00 


-5.5 


9.1 






HIP 79124 


AOV 


2 


0.6 


1.06 


0.5 


11.3 






HIP 79156 


AOV 


2 


0.5 


2.83 


6.4 


11.0 


Y 


debris 


HIP 79250 


A3III/IV 


2 


0.2 


1.29 


-19.6 


9.5 






HIP 79366 


A3V 


2 


0.6 


0.99 


2.7 


13.3 






HIP 79374 


B2IV 


2 


0.9 


0.95 


0.5 


0.9 






HIP 79392 


A2IV 


2 


0.5 


0.96 


3.0 


16.0 






HIP 79404 


B2V 


2 


0.3 


0.96 


-6.9 


3.1 






HIP 79410 


B9V 


2 


0.6 


3.77 


6.9 


10.2 


Y 


debris 


HIP 79439 


B9V 


2 


0.6 


1.99 


25.5 


12.2 


Y 


debris 


HIP 79462 


G2V 


2 


0.3 


1.51 


-3.5 


9.9 


Y 


debris 


HIP 79530 


B6IV 


2 


0.4 


0.97 


1.3 


5.3 






HIP 79606 


F6 


2 


0.9 


0.96 


-9.8 


11.7 






HIP 79643 


F2 


2 


0.6 


1.85 


1.0 


19.0 


Y 


debris 


HIP 79644 


F5 


2 


0.8 


0.93 


-51.4 


30.5 






HIP 79733 


AlmA9-F2 


1 


1.4 


1.00 


-3.6 


18.1 






HIP 79739 


B8V 


2 


1.0 


0.99 


-14.4 


18.1 






HIP 79771 


B9V 


2 


1.1 


0.97 


-19.7 


32.5 






HIP 79785 


B9V 


2 


0.2 


1.07 


-3.4 


6.9 






HIP 79860 


AOV 


1 


0.3 


1.03 


-1.3 


16.7 






HIP 79878 


AOV 


1 


0.0 


3.23 


23.1 


14.3 


Y 


debris 


HIP 79897 


B9V 


2 


0.5 


1.11 


1.3 


8.0 






HIP 80024 


B9II/III 


2 


0.7 


3.36 


-8.7 


65.4 


Y 


debris 


HIP 80059 


A7III/IV 


2 


0.5 


0.99 


-21.9 


14.1 






HIP 80088 


A9V 


2 


0.4 


3.75 


122.4 


24.6 


Y Y 


debris 


HIP 80112 


B1III 


2 


1.2 


0.95 


2.7 


6.4 






HIP 80130 


A9V 


2 


0.6 


1.04 


2.3 


12.7 






HIP 80311 


AOV 


1 


1.0 


0.86 


10.1 


50.5 






HIP 80324 


AOV+AOV 


1 


0.0 


1.03 


5.4 


11.3 






HIP 80338 


B8II 


2 


2.3 


1.42 


62.9 


17.2 






HIP 80493 


B9V 


1 


0.9 


1.01 


4.2 


18.3 






HIP 80896 


F3V 


1 


0.0 


0.97 


30.4 


12.1 







Debris Disks in Upper Sco 

TABLE 2 — Continued 
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Source 



SpT 



HIP 81266 








B0V 


HIP 82319 








F3V 


HIP 82397 








A3V 




PBB2002 


USco 


J155624 


8- 


222555 


M4 




PBB2002 


USco 


J155625 


7- 


224027 


M3 




PBB2002 


USco 


J155629 


5- 


225657 


M3 




PBB2002 


USco 


J155655 


5- 


225839 


M0 




PBB2002 


USco 


J155706 


4- 


220606 


M4 




PBB2002 


USco 


J155729 


9- 


225843 


M4 




PBB2002 


USco 


J155746 


6- 


222919 


M3 




PBB2002 


USco 


J155829 


8- 


231007 


M3 




PBB2002 


USco 


J155918 


4- 


221042 


M4 




PBB2002 


USco 


J160159 


7- 


195219 


M5 




PBB2002 


USco 


J160210 


9- 


200749 


M5 




PBB2002 


USco 


J160222 


4- 


195653 


M3 




PBB2002 


USco 


J160226 


2- 


200241 


M5 




PBB2002 


USco 


J160245 


4- 


193037 


M5 




PBB2002 


USco 


J160329 


4- 


195503 


M5 




PBB2002 


USco 


J160341 


8- 


200557 


M2 




PBB2002 


USco 


J160343 


3- 


201531 


M2 




PBB2002 


USco 


J160350 


4- 


194121 


M5 




PBB2002 


USco 


J160357 


9- 


194210 


M2 




PBB2002 


USco 


J160418 


2- 


191055 


M4 




PBB2002 


USco 


J160428 


4- 


190441 


M3 




PBB2002 


USco 


J160435 


6- 


194830 


M5 




PBB2002 


USco 


J160439 


1- 


194245 


M4 




PBB2002 


USco 


J160449 


9- 


203835 


M5 




PBB2002 


USco 


J160456 


4- 


194045 


M4 




PBB2002 


USco 


J160502 


1- 


203507 


M2 




PBB2002 


USco 


J160508 


3- 


201531 


M4 




PBB2002 


USco 


J160516 


1- 


193830 


M4 




PBB2002 


USco 


J160517 


9- 


202420 


M3 




PBB2002 


USco 


J160521 


9- 


193602 


Ml 




PBB2002 


USco 


J160525 


5- 


203539 


M5 




PBB2002 


USco 


J160528 


5- 


201037 


Ml 




PBB2002 


USco 


J160531 


3- 


192623 


M5 




PBB2002 


USco 


J160532 


1- 


193315 


M5 




PBB2002 


USco 


J160545 


4- 


202308 


M2 




PBB2002 


USco 


J160600 


6- 


195711 


M5 




PBB2002 


USco 


J160611 


9- 


193532 


M5 




PBB2002 


USco 


J160619 


3- 


192332 


M5 




PBB2002 


USco 


J160622 


8- 


201124 


M5 




PBB2002 


USco 


J160628 


7- 


200357 


M5 




PBB2002 


USco 


J160643 


8- 


190805 


K6 




PBB2002 


USco 


J160647 


5- 


202232 


M2 




PBB2002 


USco 


J160702 


1- 


201938 


M5 




PBB2002 


USco 


J160704 


7- 


201555 


M4 




PBB2002 


USco 


J160707 


7- 


192715 


M2 




PBB2002 


USco 


J160708 


7- 


192733 


M4 




PBB2002 


USco 


J160719 


7- 


202055 


M3 




PBB2002 


USco 


J160739 


4- 


191747 


M2 




PBB2002 


USco 


J160801 


4- 


202741 


K8 




PBB2002 


USco 


J160801 


5- 


192757 


M4 




PBB2002 


USco 


J160802 


4- 


202233 


M5 




PBB2002 


USco 


J160804 


3- 


194712 


M4 




PBB2002 


USco 


J160815 


3- 


203811 


M3 




PBB2002 


USco 


J160818 


4- 


190059 


M3 




PBB2002 


USco 


J160823 


2- 


193001 


K9 




PBB2002 


USco 


J160823 


5- 


191131 


M2 




PBB2002 


USco 


J160823 


8- 


193551 


Ml 




PBB2002 


USco 


J160825 


1- 


201224 


Ml 




PBB2002 


USco 


J160827 


5- 


194904 


M5 




PBB2002 


USco 


J160843 


1- 


190051 


M4 




PBB2002 


USco 


J160854 


0- 


203417 


M4 




PBB2002 


USco 


J160900 


0- 


190836 


M5 




PBB2002 


USco 


J160900 


7- 


190852 


K9 




PBB2002 


USco 


J160903 


9- 


193944 


M4 




PBB2002 


USco 


J160904 


0- 


193359 


M4 




PBB2002 


USco 


J160913 


4- 


194328 


M3 




PBB2002 


USco 


J160915 


8- 


193706 


M5 




PBB2002 


USco 


J160933 


8- 


190456 


M2 




PBB2002 


USco 


J160946 


4- 


193735 


Ml 




PBB2002 


USco 


J160953 


6- 


175446 


M3 




PBB2002 


USco 


J160954 


4- 


190654 


Ml 



Ref 
SpT 


A v 




Ratio 1 




Excess? 


TvDe 


2 4 (i m 


70^tm 


cr(70/im) 


24/im 70/^m 


1 


0.3 


1.17 


-0.3 


1.7 






2 


0.0 


1.00 


10.5 


14.2 






1 


0.0 


2.06 


-4.9 


12.8 


Y 


debris 


7 


1.7 


22.18 


110.8 


82.8 


Y 


primordial 


7 


1.0 


1.06 


23.1 


139.6 






7 


0.9 


1.06 


50.6 


123.2 






7 


0.7 


1.06 


13.7 


48.3 






7 


2.0 


10.36 


184.4 


125.3 


Y 


primordial 


7 


1.4 


11.21 


433.6 


219.2 


Y 


primordial 


7 


1.2 


1.21 


-6.2 


153.2 






7 


1.3 


28.99 


76.5 


120.7 


Y 


primordial 


7 


1.3 


1.01 


12.0 


73.1 






6 


0.6 


0.99 


231.4 


237.9 






6 


0.8 


0.84 


-224.6 


301.4 






6 


1.2 


1.05 


99.6 


132.4 






6 


0.3 


0.90 


21.7 


166.3 






6 


1.1 


0.97 


-1.9 


211.8 






6 


0.3 


1.08 


-2.0 


213.6 






6 


0.9 


1.03 


-16.6 


51.4 






6 


0.9 


0.99 


39.7 


80.6 






6 


0.6 


1.00 


-34.4 


108.8 






6 


1.7 


22.45 


219.3 


86.7 


Y 


primordial 


6 


0.8 


0.91 


-253.6 


198.3 






6 


1.0 


1.00 


-12.9 


36.9 






6 


0.8 


1.11 


10.0 


157.1 






6 


0.1 


1.48 


-33.1 


141.9 


Y 


debris 


6 


0.7 


1.07 


-82.4 


142.1 






6 


0.2 


1.03 


141.7 


134.1 






6 


1.8 


1.02 


20.1 


53.5 






6 


0.4 


1.05 


-73.6 


57.0 






6 


0.7 


0.89 


-65.6 


325.6 






6 


0.6 


0.94 


15.7 


37.8 






6 


1.2 


1.03 


56.8 


87.5 






6 


0.8 


10.22 


17.4 


176.4 


Y 


primordial 


6 


1.0 


1.06 


-35.4 


86.4 






6 


0.5 


1.05 


122.6 


240.0 






6 


0.6 


5.93 


271.7 


142.3 


Y 


primordial 


6 


2.2 


26.32 


298.4 


113.8 


Y 


primordial 


6 


0.6 


12.67 


143.4 


71.6 


Y 


primordial 


6 


1.5 


0.97 


52.2 


181.1 






6 


0.6 


0.85 


-496.6 


312.1 






6 


0.2 


37.40 


289.1 


214.0 


Y 


primordial 


6 


0.0 


1.05 


153.5 


105.6 






6 


1.8 


5.29 


46.7 


33.6 


Y 


primordial 


6 


1.4 


1.09 


-54.9 


70.1 






6 


1.0 


14.53 


112.6 


187.0 


Y 


primordial 


6 


0.8 


0.96 


-59.3 


296.9 






6 


1.8 


1.18 


-30.0 


53.2 






6 


1.2 


1.37 


113.0 


190.6 


Y 


debris 


7 


1.6 


1.44 


2.3 


164.7 


Y 


debris 


6 


1.4 


2.46 


53.8 


72.2 


Y 


debris 


6 


1.5 


1.16 


2.5 


57.2 






6 


0.3 


1.00 


-24.1 


56.4 






6 


0.7 


1.02 


20.4 


129.0 






6 


0.0 


0.84 


-200.1 


218.7 






6 


1.5 


1.22 


-206.0 


140.2 






7 


1.3 


1.09 


74.1 


88.9 






6 


1.5 


33.47 


405.3 


59.2 


Y Y 


primordial 


6 


1.1 


1.18 


49.9 


83.2 






6 


1.5 


1.01 


36.6 


35.7 






6 


1.4 


1.10 


24.4 


84.4 






6 


1.1 


10.56 


174.1 


108.8 


Y 


primordial 


6 


1.1 


1.01 


-146.1 


96.6 






7 


1.3 


1.02 


72.1 


106.4 






7 


0.7 


30.44 


81.9 


158.3 


Y 


primordial 


6 


0.8 


101.78 


834.3 


93.0 


Y Y 


primordial 


6 


0.7 


0.88 


24.4 


105.9 






7 


1.6 


1.00 


13.7 


149.9 






6 


1.4 


0.99 


-93.8 


140.3 






6 


0.3 


0.97 


184.7 


296.9 






6 


1.7 


1.08 


-142.9 


89.8 






6 


1.6 


1.07 


13.8 


76.6 






7 


1.9 


22.76 


341.5 


252.8 


Y 


primordial 


6 


0.9 


1.52 


133.9 


92.1 


Y 


debris 
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TABLE 2 — Continued 



PBB2002 
PBB2002 
PBB2002' 
PBB2002' 
PBB2002' 
PBB2002 
PBB2002 
PBB2002 
PBB2002 
PBB2002 
PBB2002 
PBB2002 
PPM 732705 
PPM 747651 
PPM 747978 



Source 

USco J160959.4- 
USco J161010.4- 
USco J161011.0- 
USco J161014.7- 
USco J161021.5- 
USco J161024.7- 
USco J161026.4- 
USco J161031.9- 
USco J161052.4- 
USco J161115.3- 
USco J161118.1- 
USco J161420.2- 



180009 
194539 
194603 
191909 
194132 
191407 
193950 
191305 
193734 
175721 
175728 
190648 



PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 



RX J1541.1- 
RX J1548.0- 
RX J1550.0- 
RX J1550.9- 
RX J1554.0- 
RX J1558.1- 
RX J1600.6- 
RX J1600.7- 
RX J1602.8- 
RX J1603.6- 
SAO 183706 
ScoPMS 13 
ScoPMS 17 
ScoPMS 21 
ScoPMS 23 
ScoPMS 27 
ScoPMS 28 
ScoPMS 29 
ScoPMS 31 
ScoPMS 32 
ScoPMS 45 
ScoPMS 52 



J153557.8- 
J154413.4- 
J155106.6- 
J155716.6- 
J155750.0- 
J155812.7- 
J155847.8- 
J160000.7- 
J160013.3- 
J160031.3- 
J160042.8- 
J160108.0- 
J160147.4- 
J160200.3- 
J160239.1- 
J160251.2- 
J160354.9- 
J160357.6- 
J160421.7- 
J160539.1- 
J160612.5- 
J160654.4- 
J160703.9- 
J160814.7- 
J160831.4- 
J160856.7- 
J160930.3- 
J161302.7- 
J161318.6- 
J161329.3- 
J161402.1- 
J161411.0- 
J161459.2- 
J161618.0- 
J161933.9- 
2656 



232405 
252258 
240218 
252918 
230508 
232835 
175800 
250941 
241810 
202705 
212737 
211318 
204945 
222123 
254208 
240156 
203137 
203105 
213028 
215230 
203647 
241610 
191132 
190833 
180241 
203346 
210459 
225744 
221248 
231106 
230101 
230536 
275023 
233947 
222828 



2908 

2312 

2534 

2920 

2405 A 

2159 

2343 

2401A 

2245 



SpT 

M4 

M3 

M5 

M3 

M3 

M3 

M4 

K7 

M3 

Ml 

M4 

K5 

G6 

G3 

G3 

K3: 

Ml 

M2 

M0 

M0 

G2 

K3 

GO 

M0 

Ml 

K7 

M0 

M0 

Ml 

K7 

K4 

M0 

K5 

K2 

Ml 

K5 

M3 

Ml 

K2 

M0 

K5 

M0 

K4 

G9 

Kl 

G4 

K0 

G5 

G7 

K0 

G7 

G9 

M2 

F9 

M0 

K4 

G9 

M2 

K0 

G9 

G8e 

M1.5V 

M1V 

K1IV 

K5IV 

K2IV 

M1V 

M2V 

M0.5V 

M3V 

K5IV 

K0IV 



Ref 
SpT 


Ay 




Ratio 1 




Excess? 


Type 


24/im 


70/im 


cr(70/im) 


24/im 70/im 


6 


0.7 


96.14 


1006.3 


156.3 


Y Y 


primordial 


6 


1.0 


1.09 


74.4 


109.7 






7 


0.5 


2.02 


58.2 


302.8 


Y 


debris 


6 


0.9 


1.50 


24.0 


91.5 


Y 


debris 


7 


1.6 


1.12 


-96.2 


68.7 






6 


1.5 


0.91 


86.3 


145.4 






7 


1.9 


0.94 


-38.2 


111.6 






6 


1.1 


1.04 


16.3 


62.2 






7 


2.3 


1.96 


78.0 


92.8 






6 


1.6 


17.71 


41.1 


27.1 


Y 


primordial 


6 


0.9 


0.99 


-61.3 


48.0 






7 


1.8 


40.07 


201.1 


21.6 


Y Y 


primordial 


4 


0.7 


1.00 


-17.3 


29.2 






4 


0.3 


1.05 


36.2 


22.7 






4 


0.7 


0.94 


3.0 


50.6 






5 


0.7 


0.90 


5.5 


66.9 






5 


0.6 


1.00 


-49.0 


52.9 






3 


0.4 


1.16 


-103.2 


73.5 






3 


0.0 


0.94 


-12.6 


48.9 






3 


0.1 


0.95 


36.0 


50.6 






3 


0.4 


2.72 


54.7 


75.0 


Y 


debris 


5 


1.8 


1.41 


-46.9 


32.1 


Y 


debris 


3 


0.0 


1.09 


47.9 


55.1 






3 


0.6 


1.50 


-74.1 


58.8 


Y 


debris 


3 


1.0 


0.98 


-5.8 


30.5 






3 


0.8 


1.01 


6.5 


37.5 






3 


0.0 


0.99 


-0.9 


26.6 






3 


0.8 


1.01 


-7.8 


26.6 






3 


1.1 


1.48 


-31.2 


32.5 


Y 


debris 


5 


0.0 


0.96 


20.9 


72.2 






3 


0.6 


1.22 


-41.0 


47.1 






3 


0.9 


0.99 


-1.5 


30.6 






3 


0.9 


33.16 


176.9 


21.3 


Y Y 


primordial 


3 


1.0 


58.15 


6204.3 


623.8 


Y Y 


primordial 


3 


0.8 


1.17 


-10.4 


61.7 






3 


1.8 


0.98 


1.6 


35.3 






5 


0.0 


0.93 


3.0 


33.4 






3 


1.1 


1.05 


-23.7 


36.4 






5 


1.6 


0.97 


15.9 


32.9 






5 


0.1 


1.06 


28.1 


47.7 






5 


1.4 


0.97 


-32.7 


29.7 






5 


0.0 


1.14 


59.3 


40.2 






5 


2.3 


1.00 


-17.9 


24.7 






5 


0.9 


0.98 


-20.7 


12.2 






5 


2.0 


1.03 


28.8 


25.4 






5 


2.0 


1.08 


-50.0 


33.4 






5 


2.4 


18.02 


49.1 


8.0 


Y Y 


primordial 


5 


1.0 


1.57 


85.8 


81.1 


Y 


debris 


5 


1.1 


1.11 


-25.8 


24.6 






5 


1.3 


1.03 


-22.3 


27.5 






5 


1.6 


1.03 


-152.6 


56.6 






3 


0.7 


1.05 


-12.7 


39.8 






3 


0.7 


0.97 


3.0 


56.5 






3 


0.2 


0.91 


8.8 


34.7 






3 


0.4 


0.95 


10.4 


27.9 






3 


0.5 


0.90 


-190.8 


156.1 






5 


0.5 


1.11 


21.1 


33.6 






3 


0.5 


0.90 


-94.6 


129.1 






3 


2.3 


0.99 


2.1 


24.3 






3 


1.2 


1.00 


9.8 


29.1 






4 


1.7 


1.04 


-15.3 


21.5 






8 


0.2 


0.95 


-13.7 


44.9 






8 


0.7 


14.68 


12.7 


3.3 






5 


0.2 


0.98 


-21.0 


67.3 






8 


0.2 


0.92 


48.9 


23.4 






5 


0.2 


0.94 


12.6 


19.7 






8 


0.7 


0.96 


11.6 


68.5 






8 


0.7 


1.02 


-10.5 


36.2 






8 


0.9 


48.89 


404.1 


43.2 


Y Y 


primordial 


8 


0.6 


0.92 


-10.9 


106.6 






8 


0.7 


1.44 


71.2 


51.7 


Y 


debris 


8 


1.3 


0.97 


35.3 


24.6 
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Source 


SpT 


Ref A v 


Ratio 1 


Excess? 


Type 






SpT 24>m 


7U^m o-(7U^m) 


24/im 7U/^m 





References. (1) IHouU dl982l): f2 ) I Houk fc Smith-Moorel dl988l): ( 3) IKunkell <199Sfl: (4) E . Mamajek, private communication; (5) 

IPreibisch et al.l 1 119981) ; (6) IPreibisch et al.lf2001l) ; ( 7) IPreibisch et al.l 1 120021) ; (8) IWalter et al. I 1 11 9941) 

^ Ratio of the observed flux density to the photospheric value. The 24/im ratios have an adopted uncertainty of 0.07. 
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TABLE 3 
Disk Types in Upper Sco 



Spt 


N (stars) 




N(disks) 


Be 


Debris .Primordial 


B 


37 


2 


6 


A 


25 





7 


F 


9 





1 


G 


21 





3 


K 


27 





2 7 


M 


101 





9 17 


Total 


220 


2 


28 24 
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TABLE 4 

Compilation of Debris Disk Surveys 



Region 


Aee 


Aee 


Nstars 


Minimum excess 3. 




SpT 


Photometry 




(NLyv) 


Reference 




(percent J 


ar y 


D Limit c 


.R cfcrcncc 


CrA 


2-16 


21 


7 


10 


G6 


K5 


4 


IC 348 


2 


17 


6 


30 


AO 


Al 


15 


a Ori 


2.5 


29 


11 


39 


B2 


A3 


13 


Upper Sco 


5 ± 1 


25 


220 


32 


B0 


M5 


34 






5 


30 


F0 


G5 


5 








16 


10 


GO 


K3 


4 








3 


8 


B9 


B8 


32 


Orion OB lb 


5 ± 1 


3 


21 


31 


Bl 


A8 


12 


r) Cha 


6 ± 2 


18 


16 


31 


B8 


M5.75 


9 






1 


8 


B8 


B8 


5 


TW Hya 


8 


37 


21 


60 


AO 


M5 


16 


Orion OBla 


8.5 ± 1.5 


3 


21 


32 


Bl 


F0 


12 


P Pic 


12 ± 6 


36 


30 


27 


AO 


M4.5 


27 








3 


10 


F8 


K7 


4 








2 


10 


AO 


A5 


32 


UCL 


16 


26 


11 


30 


F0 


G8 


5 








23 


10 


F8 


K3 


4 








16 


8 


B7 


A2 


32 


LCC 


17 


26 


24 


30 


F0 


G5 


5 








11 


10 


GO 


K2 


4 


NGC 2232 


29 ± 4 


19 


15 


30 


Bl 


A8 


7 


Tucana-Horologium 


30 


20 


9 


27 


B8 


M3 


27 


NGC 2547 


38 ± 5 


22 


31 


15 


B9 


F5 


11 


IC 2391 


50 ± 5 


1 


29 


15 


B3 


K4 


30 


IC 2602 d 


50 ± 5 


1,31 


7 


8 


B6 


B9 


32 








5 


10 


GO 


K0 


4 


NGC 2422 


80 ± 20 


35 


31 


15 


B5 


A 


28 


a Per 


85 ± 10 


1 


13 


10 


F5 


K2 


4 








5 


8 


B8 


AO 


32 


Pleiades 


130 ± 20 


1 


20 


10 


F6 


K2 


4 








72 


15 


B6 


K3 


10 








6 


8 


B8 


AO 


32 


NGC 2516 


125 ± 25 


19 


51 


15 


B8.5 


F0 


28 


Ursa Major 


500 ± 100 


11 


7 


8 


AO 


A6 


32 


Coma Berenices 


500 ± 100 


23 


5 


8 


A2 


Al 


32 


Hyades 




24 


7Q 
/ 


1 8 
±0 


A 1 
ill 


k'7 

IV t 


r 
u 








22 


10 


F8 


K2 


4 








11 


8 


Al 


A7 


32 


Praesepe 


757 ± 36 


8 


90 


15 


Al 


G6 


8 








5 


8 


AO 


A6 


32 


Field 


150-4570 


2 


38 


30 


F0 


K5 


2 




3-4000 


21 


186 


10 


F5 


K7 


4 




5-680 


28 


30 


15 


B7 


A5 


28 




5-850 


32 


85 


8 


B7 


A7 


32 




190-5000 


33 


98 


10 


F0 


K4.5 


33 



References. — (1) IBarrado v Navascues et all J20O4D: ( 2) Boichman ct al. (2 0061): (3) IBriceno et al.l 120051); ( 4) [Carpenter et al. (2008): 

(5) Chen ct al. (2005): (6) Cicza ct al. (2008); (7) Curric ct al. (2008b): (8) Caspar ct al. (2009); (9) Gauticr ct al. (2008): (10) Gorlova ct al. 

(2006) ; (11) IGorlova et al.l 120071); f 12) IHernan.dez et al.l l|2006D; (1 31 Hernan dez et al.l (120071); (14)IKing et a l. (2003); ( 15) ILada et al.l (120061) ; 

(16)ILow et alj (120051) ; (1 7)ILuhman et al l (120031) ; (18) [Luhman &; Steeghsl fl2004) ; (19) ILvra et alj (I20qd); (20 ) IMamaiek et al.l (120041); (21 ) 
IMever et al.l 1120061)7 (22) [ Navlor &: Jeffried H2006I): (23) (Odenkirchen et al.l U998I); (24) (Perrvman et all 111993) ; (2 5) IPreibisch et al. | 120021) : 

(26) IPreibisch fc MamaieM <2008l): (27) fR ebull et alt 1)20081); (28) IRieke et al J 1 12003) ; (29) ISherrv et aU 11200811: (301 ISiegler et al.l (120071); "(31) 
IStauffer et alj U997T); (32)ISu et alj dgOOfii) : (33) ITrilling ct al] ( 120081) : (34) This study; (35)|van RcnsbciEcn ct aLI dl978l) ; (36) IZuckernian et all 

(120011) ; (37) IZuckerman fc Sond (120041) 

a The minimum percent 24/im flux density above the photosphere used to identify an excess 
k Earliest spectral type in the survey 

c Latest spectral type where the minimum excess could be detected 

We adopt an age of 50 Myr for IC 2391, and assume IC 2391 and IC 2602 have the same age following Stauffcr et alj 1)19971) . 
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( J - K s)o 

Fig. 1. — 24^tm to 4.5/im flux density ratio (-R24/4.5) as a function of dereddened J — K s color. The solid line indicates the best-fit 
line [log -R24/4.5 = 0.098 (J — K B ) — 1.40] to the main locus of points that is assumed to represent the stellar photosphere. The dashed 
lines indicate the ±4<r limits about the fit (<r = 0.03 dex). Black circles indicate sources with an 8fim or Wfim excess from Paper I. Black 
stars are sources with detectable excesses at wavelengths > 2Afim, while gray circles indicate sources where the colors are consistent with 
photospheric emission. Sources with questionable photometry from nebulosity or a contaminating source are indicated by open circles, and 
known Be stars are marked by the crosses. The source [PBB2002] USco J161420. 2— 190648 is offscale on this plot at (J — A" s ) ,-R24/4.5 = 
2.14, 0.27. 
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HIP 77859 





HIP 78207 




HIP 7791 1 




[PBB2002] USco J1 55624.8-222555 





Fig. 2. — 24/im images for the 54 Upper Sco stars with a 24/im excesses. The left panel for each star shows a 1' X 1' region centered on 
the stellar position. The right panel shows the PRF-subtracted image, where all identified stars in the field of view have been fitted and 
removed. All images are displayed in a logarithmic scale. 
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[PBB20021 USco J 155706.4-220606 [PBB20021 USco J 155729.9-225843 






[PBB2002] USco J1 55829.8-231 007 [PBB2002] USco J1 60357.9- 1 9421 




• 




[PBB2002] USco J1 60439.1 -1 94245 [PBB2002] USco J1 60525.5-203539 







[PBB2002] USco J1 60532.1 -1 9331 5 [PBB2002] USco J 1 60545.4-202308 







[PBB2002] USco J1 60600.6-1 9571 1 [PBB2002] USco J1 60622.8-201 1 24 







[PBB2002] USco J 1 60643.8-1 90805 [PBB2002] USco J1 60702.1 -201 938 





1002] USco J 




[PBB2002] USco J 1 60708.7-1 92733 [PBB2002] USco J1 6071 9.7-202055 







[PBB2002] USco J1 60739.4-1 91 747 [PBB20021 USco J1 60823.2-193001 






[PBB2002] USco J 1 60827.5-1 94904 [PBB2002] USco J 1 60900.0-190836 







Debris Disks in Upper Sco 
[PBB2002] USco J1 60900.7- 1 90852 [PBB2002lUScoJ1 60953.6-1 75446 
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* 



160900.7-190852 [PBB2002] USco J 

■ ■ 




[PBB2002] USco J1 60954.4- 1 90654 [PBB2002] USco J1 60959.4-1 80009 





[PBB2002] USco J1 61 01 1 .0- 1 94603 [PBB2002] USco J1 61 01 4.7-1 9 1 909 







[PBB2002] USco J1 61 1 1 5.3-1 75721 [PBB2002] USco J1 61 420.2-1 90648 




[PZ99] J1 5581 2.7-232835 





[PZ99] J160013. 3-241810 





[PZ99] J1 60200.3-2221 23 
[PZ99] J1 60421. 7-21 3028 



* 




[PZ99] J160357. 6-203105 



* 




ScoPMS 31 



[PZ99] J1 61 41 1.0-230536 



ScoPMS 45 





2G 



50 - 



40 - 



30 - 
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20 - 



10 - 
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[S I/ (70 M m) observed - S„(70 M m) photosphere ] / a 



Fig. 3. — Histogram of the signal-to-noise ratio for the observed 70/xm excess above the stel lar p hotosphere. The histogram excludes 
two stars (ScoPMS 17 and HIP 80338) for which the 70/^m photometry is suspect (see Section |2,2[ I, Sources were identified as having a 
significant excess if the observed 70/xm flux density exceeded the photospheric flux density by > 3<x. 
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A v = 20 
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Fig. 4. — The dereddened 2Afim to 4.5/mi flux ratio as a function of the dereddened 8fim to 4.5/xm flux ratio for B7-A9 stars (top), 
F0-G9 stars (middl e), and K0-M5 sta rs (bottom). Filled circles represent stars in Upper Sco, and crosses indicate stars in IC 348 that have 
a 24/xm detection (Lada ct al. 2006). Black symbols indicate stars with an infrared excess at 24/jm (excluding two Be stars), and gray 
symbols indicate stars that lack a 24/jm excess. Crosses bounded by circles and squares are stars in IC 348 that have been classified by 
ILada et all (2006) as "anemic" or "flat" based on the slope of the SED in the Spitzer IRAC bands. The dashed lines indicate the range of 
slopes (a = —2.56 to 1.80) used to define anemic disks (a = dlog(vS v )/dlog(v)). A few anemic disks lie outside these boundaries since the 
actual slopes were computed using all 4 IRAC bands. A GO star in iC 348 is offscale at (6.4, 4.8). 
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Fig. 5. — The 16^tm (top panel) and 24^tm (bottom panel) excess in Upper Sco normalized by the photospheric emission as a function 
of the dereddcncd J — K B color. Symbols indicate debris (filled circles) and primordial (asterisks) disks. Triangles indicate the 3cr upper 
limits for stars without detected excesses. Be stars and sources with contaminated photometry have been omitted. The solid curves show 
the emission from model debris disks with a power-law distribution of particle sizes at a given orbital radius. The models were normalized 
to the debris properties around B/A stars, and computed around stars of various temperatures after factoring differences in stellar heating 
and radiation blowout size. Open circles indicate the nine B/A stars used to normalize the models. See text for explanation of model 
behavior as a function of photospheric color. The results indicate that the magnitude of the 16/im and 24/im excess around the debris disks 
can be explained by planetesimals belts that have a similar range of orbital radii across all spectral types. 
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Fig. 6. — Fraction of sources with a 24/im excess from a debris disk versus stellar age for the regions summarized in Table [4] The plots 
in the left panels show the results for a 24/xm excess detection threshold of 15%, and the right panels for 32%. Results are presented for 
three spectral type ranges. Grey circles represent individual clusters or associations that contain at least 5 stars in the appropriate mass 
range. Black circles represent the ensemble average of field stars, clusters, and associations regardless of the sample size. 
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Fig. 7. — Ratio of the 24/jm excess normalized by the photospheric 24^tm flux density versus stellar age compiled from the regions 
summarized in Table [4] Filled circles represent stars with a detected 24/im excess that are likely debris disks, and open triangles represent 
3<t upper limits for stars without excesses. Asterisks indicate optically thick disks (including "anemic" and "evolved" disks). Known Be 
stars and sources with c ontaminated 24/xm photometry have been omitted. The dashed line shows the upper envelope to the 24^tm excess 
fromfRickc ct al. (2003) for ages older than 5 Myr, and extrapolated here to younger ages. 
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Fig. 8. — Expected temporal variation in the 24/^m excess relative to the stellar photosphere for AO (solid black curve), F0 (dashed), 
GO (dotted), K0 (dash-dot), and M0 stars (dash-dot-dot-dot) stars. The 24/xm excess for a given curve has been normalized by the peak 
excess. The stellar mass and luminosity will vary with age for a fixed spectr al type as the star evolves toward the main-sequence. The 
stellar parameters for a given age and spectral type were estimated from the Sicss et al. (2000) pre-main-sequence models. The model 
planctesimal belt is at an orbital radius of 15 AU. The debris emission varies in time as stars evolve, which changes the stellar heating 
and the radiation blowout size. These results show that if the mass of the planctesimal belt is fixed, A-K stars younger than 10 Myr will 
have a lower 24/xm excess since the smaller grains have been blown out of the system by the greater stellar luminosity. For comparison, 
the solid gray curve shows the variation in the normalized 24/^m excess from a planetesimal belt with a mass of the minim um mass solar 
nebula that extends between 30 and 150 AU around a 2Mg star with constant stellar luminosity (Kcnvon & Bromlcv 2008). 
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Fig , 9. — The 24/j m excess for M0-M5 star s as a function of stellar age. Exce ss sources are shown for IC 348 (age 2 Myr; Lada et al.l 
2006h. a Ori (2 .5 Myr: IHernandez et al1l2007l. Upper Sco ( 5 Myr), jj Cha (6 Mvr: IGautier et al.ll200l) . TW Hydra moving group (8 Myr; 
Low et alj|2005ft , NGC 2547 (38 Myr; IForbrich et alj|200a) . and IC 2391 (50 Mvr; ISiegler et alJl200W) . Asterisks indicate stars where the 
24/tm excess is thought to originate from an optically thick disk; circled asterisks are optically thick disks further subdivided as "anemic" or 
"evolved" disks. Filled circles indicate stars classified as debris systems. M-stars with upper limits at 24^m have been omitted for clarity. 
Not all surveys shown here are sensitive to the stellar photosphere, and this diagram should not be used to infer the evolution of debris 
emission around M-type stars. 
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